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ABSTRACT 


The  time-harmonic  analysis  of  three  boundary  value  problems  con¬ 
taining  semi-infinite  boundaries  is  presented.  The  first  problem  considered 
is  a  parallel  plate  waveguide  with  one  plate  truncated  and  radiating  into  free 
space.  The  excitation  of  a  dielectric  slab  and  the  excitation  of  an  isotropic, 
incompressible,  plasma  slab  by  means  of  a  parallel  plate  waveguide  v/ith 
one  plate  truncated  are  the  second  and  third  problems  analyzed,  respectively. 
Both  TE  and  TEM  polarizations  are  considered  in  these  open-region  problems. 

A  function  of  a  complex  variable  is  factored  in  each  of  these  Wiener- 
Hopf  type  boundary  value  problems.  The  function  is  analytic  in  a  strip  and  is 
factored  into  a  product  of  two  functions.  One  of  these  functions  is  analytic  in 
a  half-plane  while  the  other  is  analytic  in  the  adjacent  half-plane  with  an  over¬ 
lap  in  the  regions  of  analyticity  coinciding  with  the  strip.  This  factorization  is 
obtained  by  a  technique  developed  in  this  work. 

The  technique  obtains  the  factorization  for  the  open- region  problem 
from  a  function  and  its  factorization  that  occurs  in  a  related  closed- region 
problem.  A  closed-region  problem  is  one  whose  transverse  dimensions  are 
finite.  The  chosen  closed-region  boundary  value  problem  yields  a  function  of 
a  complex  variable  which  can  be  factored.  The  factorization  of  the  function 
for  the  open-region  boundary  value  problem  is  obtained  by  taking  the  limit,  as 
a  parameter  approaches  infinity,  of  the  function  and  factorization  appropriate 
to  the  closed- region  structure.  By  this  means  the  factorization  and  hence  the 
solution  to  the  open- region  boundary  value  problem  is  obtained. 


It  is  also  founu  mat  the  limiting  procedure  may  be  used  to  obtain  more 
than  just  the  open- region  factorization.  It  is  shown  that  the  limit  of  the  com¬ 
plete  closed-region  solution  becomes  the  open-region  solution.  Hence,  this 
yields  one  possible  method  for  the  solution  of  problems  of  this  type. 

The  results  of  the  numerical  computations  are  presented.  These  in¬ 
clude  the  average  power  reflected  in  the  waveguide,  the  average  power  radiated 
in  the  space  wave,  the  average  power  transmitted  by  the  surface  waves,  and 
the  radiation  pattern  of  the  space  wave. 
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1.  INTRODUCTION  }i 

l  i  \ 

"  }  t 

The  Mmo-harmoruc  analys  s  of  certain  radiation  and  diffraction  prob¬ 
lem;;  with  semi- infinite  boundaries  requires  solutions  of  the  steady- state  wave 

-s  f  \  - 

|  ! 

equation  satisfying  various  boundary  conditions.  This  class  of  problems  is  con- 

i  - 

ventionally  formulated  in  terms  of  the  Wiener-Hopf  technique,  that  is,  at  some 

I 

point  in  the  analysis  a  complex  variable  equation  is  solved  by  analytic  continua¬ 
tion,  ‘  An  exhaustive  discussion  and  numerous  illustrations  of  this  technique  may 

i 

be  found  in  Noble  [1958}. 

j 

Difficulty  v.-.\th  tb«»  Wmner-Hopf  technique  is  encountered  because  a  fac¬ 
torization  of  a  function  of  a  complex  variable  must,  be  made.  This  function  of  a 

l 

complex  variable,  hich  is  analytic  in  a  strip,  must  be  factored  into  a  product 
of  two  functions.  One  function  of  '>e  p.oduct  is  analytic  in  a  half-plane  while 

i 

the  other  is  analytic  in  'he  adjacent  half -plane,  with  an  overlap  in  the  regions  of 

j 

analyticity  coinciding  with  the  strip. 

j 

In  this  work  a  C-R  boundary  value  problem  will  refer  to  a  closed- region 
boundary  value  problem  (one  whose  transverse  dimensions  are  finite;  see  Fig.  2, 
for  example).  An  O  R  boundary  value  problem  means  an  open  region  boundary 
value  problem  (one  in  which  radiation  may  occur;  see  Fig.  1,  for  example). 

The  factorizat.on  in  the  case  of  a  C-R  boundary  value  problem  may  be 
obtained  by  U3ing  the  infinite  product  expansion  of  an  integral  (entire)  function  , 
set:  for  example  Titchmarsh  [1932],  T*us  results  from  the  fact  that  the  function 
is  a  ratio  of  two  integral  functions.  Tve  function  to  be  factored  in  the  case  of  an 


<t  Availabls  Copy 
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O-R  boundary  value  problem  contains  branch  line  singularities  and  recourse  to  a 
formal  factorization  procedure,  for  example  Noble  [1958],  may  be  made.  How¬ 
ever,  as  usually  occurs  with  such  procedures,  specific  results  are  difficult  to 
obtain  except  In  a  few  simple  cases. 

The  possibility  of  attacking  an  0-R  boundary  value  problem  through  a 
related  C-R  boundary  value  problem,  which  by  its  nature  is  easier  to  analyze, 
has  been  suggested  and  attempted  to  a  limited  extent  by  various  authors.  Noble 
[19  58 J  expressed  interest  in  knowing  how  far  the  results  for  a  parallel  plate 
duct  (semi- infinite  parallel  plate  waveguide)  enclosed  in  a  larger  parallel  plate 
waveguide,  with  a  finite  (but  large)  spacing  between  the  plates,  could  be  used  to 
approximate  the  results  near  the  mouth  of  the  parallel  plate  duct  when  radiating 
in  unbounded  space.  Talanov  [1959],  desiring  the  analysis  of  surfave  wave 
launching  in  a  dielectric  slab  backed  by  a  perfect  conductor  by  means  of  a  semi¬ 
infinite  parallel  plate  waveguide,  enclosed  the  O-R  structure  in  a  larger  parallel 
plate  waveguide.  This  procedure  reduced  the  O-R  structure  to  a  C-R  structure. 
He  then  analyzed  this  C-R  structure  and  calculated  the  desired  field  quantities 
for  increasing  values  of  the  spacing  between  the  plates  of  the  parallel  plate  wave¬ 
guide.  He  suggested  that  the  results  obtained  for  the  C-R  structure  are  in  the 
limit,  as  the  spacing  between  the  plates  of  the  parallel  plate  waveguide  becomes 
large,  the  results  of  the  O-R  problem.  Mittra  and  Karjala  [1964]  showed  that 
the  expression  for  the  reflection  coefficient  of  a  parallel  plate  duct  enclosed  in  a 
larger  parallel  plate  waveguide  yields,  in  the  limit  of  the  waveguide  walls  ap- 
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proaching  infinity,  the  expression  ior  the  reflection  coefficient  jn  the  duct  when 
radiating  into  free  space.  Mi’tra  and  vanBlaruum  [1965]  numerically  calcu¬ 
lated  the  leflection  coefficient  in  the  duct  enclosed  in  the  larger  parallel  plate 
wavegu.  de  and  showed  tha*  ’he  numerical  values  approached,  as  the  spacing  of 
the  plates  of  the  waveguide  became  large,  the  known  numerical  value  of  the  re¬ 
flection  coefficient  for  the  duct  radiating  into  free  space.  Mittra  and  Ba’es 
[1965]  used  a  limiting  procedure  to  obtain  an  extension  of  'he  function  -  theoretic 
technique  introduced  by  Whitehead  [1951],  The  limit,  as  a  dimension  became 
infinitely  large,  of  a  certain  function  that  occurs  m  a  related  C-R  problem  gave 
the  desired  unknown  function  necessary  in  ’be  O-R  problem.  The  mod*  match¬ 
ing  technique  was  used  in  that  analysis. 

The  extension  of  a  C-R  boundary  valve  problem  solution  to  yield  the 
solution  of  an  O-R  boundary  value  problem  ;s  expected  one  >,ak->s  in’o  account 
the  physical  phenomenon  occurring,  "For  example,  consider  a  source  in  a  par¬ 
allel  plate  wayeguide  where  the  medium  has  a  sligh  loss.  At.  any  location  A 
within  the  waveguide  the  field  is  mad**  up  of  two  components;  a  direct  wave  from 
the  source  and  reflected  waves  from  the  boundary.  The  magnitude  of  th»  re 
fleeted  waves  at  A,  as  the  spacing  of  the  waveguide  walls  approaches  mfin :  ’y , 
would  approach  zero  due  to  the  loss  in  ’he  medium.  Therefore,  point  A  would 
see  only  the  incident  field  in  the  limit.  That  is,  we  arc  T.  t  with  .  source  radiat¬ 
ing  in  an  unbounded  region.  The  idea  of  a  slight  loss  in  the  medium  is  not  re¬ 
strictive,  When  th?  analysis  is  completed  the  loss  is  permitted  to  be  as  email 


4 


as  desired,  in  fact,  zero.  The  inclusion  of  the  loss  is  usually  used  regardless 
of  the  method  of  solution  of  these  problems. 

This  work  determines  a  method  by  which  the  O-R  factorization  may  be 
obtained  from  a  related  C-R  problem.  The  method  involves  a  limit,  as  a  pa¬ 
rameter  approaches  infinity,  as  suggested  by  the  physics  of  these  problems. 

The  chosen  C-R  boundary  value  problem  yields  a  function  of  a  complex  variable 
which  con  be  factored.  The  factorization  of  the  function  for  the  O-R  boundary 
value  problem  is  obtained  by  taking  the  limit,  as  the  transverse  dimension  ap¬ 
proaches  infinity,  of  the  function  and  factorization  appropriate  to  the  chosen  C-R 
structure.  By  this  means  the  factorization  and  hence  the  solution  to  the  O-R 
boundary  value  problem  is  obtained. 

It  is  found  that  the  limiting  procedure  may  be  used  to  obtain  more  than 
just  the  O-R  factorization.  It  is  shown  that  the  limit  of  the  complete  C-R  solu¬ 
tion  becomes  the  O-R  solution,  Herce,  this  yields  one  possible  method  for  the 
solution  of  problems  of  this  type.  This  is  a  rather  useful  method  as  the  C-R  so¬ 
lution  is  usually  readily  obtained. 

Obviously  there  is  more  than  one  possibility  for  the  choice  of  a  C-R 
structure.  However,  the  results  obtained  for  the  O-R  structure  are  unique  since 
the  O-R  solution  is  a  limit  point  of  the  C-R  solutions.  This  result  is  expected 
from  the  physics  of  the  problem  which  implies  that  the  field  reflected  from  a 
boundary  that  is  receding  to  infinity  in  a  lossy  medium  will  be  zero  in  the  vicin¬ 
ity  of  the  source.  Hence,  the  boundary  condition  satisfied  by  the  boundary  that 
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that  recedes  to  infinity  in  the  C-R  structure  is  immaterial. 

The  firs'-  problem  discussed  is  the  analysis  of  the  fields  associated 
with  a  parallel  plate  waveguide  having  the  top  plate  terminated  (semi- infinite). 
The  solution  is  obtained  in  closed  form  and  thus  the  method  is  clearly  demon¬ 
strated.  The  factorization  is  verified  by  reference  to  the  solct.on  for  the 
fields  of  a  parallel  plate  duct  obtained  by  Noble  [1958],  since  the  function  'o  be 
factored  s  the  same  in  each  problem.  The  launching  of  surface  waves  on  a  di¬ 
electric  slab  with  a  relative  dielectric  constant  greater  than  one  by  means  of  a 
semi- infinite  parallel  plate  waveguide  is  analyzed.  Solutions  for  both  TE  and 
TEM  excitations  are  obtained.  Numerical  results  for  the  power  reflected  m 
the  waveguide  power  Tapped  in  the  surface  waves,  power  radiated  by  the  space 
wave,  and  also  the  radiation  pattern  of  the  space  wave  are  obtained  for  various 
parameters.  The  power  results  for  the  TEM  excitation  are  compared  with 
those  obt?vned  by  Angulo  and  Chang  [  1 9 c 9 1  who  worked  with  the  formal  factori¬ 
zation  procedure  and  the  differences  in  *he  results  are  noted.  The  final  prob¬ 
lem  analyzed  is  the  case  where  tt\e  dielectric  slab  is  replaced  by  an  incom¬ 
pressible,  isotropic,  plasma  slab.  This  gives  the  possibility  of  a  relative  di¬ 
electric  constant  less  than  one.  Again,  numerical  results  for  th»  power  reflec¬ 
ted  in  the  waveguide,  power  radiated  in  the  space  wave,  and  the  radiation  pat¬ 
tern  of  the  space  wave  are  presented  for  various  parameters.  No  trapped 
waves  can  occur  m  this  case. 

The  boundary  value  problems  investigated  here  are  formula'ed  by  a 
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method  used  by  Jones  [  1 9 50]  as  opposed  to  an  integral  equation  approach.  Four¬ 
ier  transforms  are  applied  directly  to  the  partial  differential  equation  and  the 
complex  variable  equation  is  obtained  without  the  use  of  an  integral  equation. 

The  integral  equation  approach  would  lead  to  equivalent  results  as  the  integral 
equation  would  be  of  the  Wiener- Hopf  type;  see  for  example  Morse  and  Fesh- 
bach  [1953].  Jones'  method  also  has  the  advantage  that  the  application  of  the 
edge  condition,  Meixner  [1954],  which  is  necessary  in  this  type  of  problem,  may 
be  clearly  applied. 

The  usual  method  of  calculating  the  far  field  is  by  means  of  saddle  point 
integration.  However,  the  problems  considered  here  are  such  that  the  far  field 
pattern  may  be  obtained  more  directly  by  using  an  equivalent  Huygen  source  in 
the  aperture.  The  far  field  pattern  is  then  related  to  the  Fourier  transform  of 
this  aperture  distribution  and  in  these  problems  becomes  an  evaluation  of  a 


function  on  an  interval. 
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2.  RADIATION  FROM  A  TRUNCATED  PARALLEL  PLATE  V'AVEGUIDE 


2.  1  Formulation  of  the  Problem 

The  first  problem  considered  is  a  parallel  plate  waveguide  with  one 
plate  truncated  {semi- infinite)  and  radiating  into  free  space.  The  structure  is 
shown  in  Fig.  1.  The  incident  field,  in  the  parallel  plate  waveguide  section  of 
the  structure,  is  taken  to  be  the  TEM  mode  with  the  magnetic  field  intensity 
parallel  to  the  walls  of  the  structure.  The  case  of  a  TE  incident  field  is  not 
discussed  in  detail  as  the  function  to  be  factored  turns  out  to  be  the  same  as  in 
the  TEM  case.  However,  the  TE  case  can  be  obtained  from  section  3.  1  by  set¬ 
ting  the  relative  dielectric  constant  (/C)  to  one. 

The  incident  field  is  therefore  the  lowest  order  TM  mode. 


*  * 


Hr  e 


(2.  1) 


We  wish  to  find  the  electromagnetic  fields  which  satisfy  Maxwell's  equation  and 
the  necessary  boundary  conditions  pertaining  to  this  structure  with  a  source  as 
given  by  (2.  1). 


Maxwell's  equations  for  a  medium  with  loss  and  the  time  convention  cho 


sen  are 

V  x  H  =  (-iou£c  +<r,)  E.  }  VxE=  i  wp.  H  (2.2] 

The  loss  is  due  to  the  conductivity  (  O',  )  of  the  medium.  From  (2.  2)  it  can  be 


e 


time  convention 


9 


shown  that  the  magnetic  field  intensity  must  satisfy 


z  — 


2  — 


H  +  k  H  *  O 


(2.3) 


where 


'/2 


k-  (u)Z^o€0  +iup0  o', )  *  k.+iki 


(2.4) 


with  k  >  0  and  k  >  0.  Since  the  incident  field  is  independent  of  the  y-coordi- 
1  Z 

nate  and  the  entire  structure  is  uniform  with  respect  to  the  y-directinn^the  total 
field  will  also  be  independent  of  y.  Therefore  the  solution  of  (Z.  3)  is  equivalent 
to  solving  the  two-dimensional  wave  equation  for  the  scalar  potential  (|>  . 


Li*  +  li*  +  (,‘i  =  o 

j  +  j  i*  +  f  * 

All  the  field  quantities  are  derivable  from  ♦  ^  by  letting 

H-  =  <f>t 


(2.  5) 


(2.6) 


E,  = 


i 

(ioj€0  -  <f,)  z 


(2.  7) 


E  y  =  - 


i 


(2.8) 


Let 


<K  =  <Pi  +  4 

where  (j> .  is  the  incident  field  and  ^  is  the  scattered  field.  Obviously 


(2.  9) 
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-itz 

*  o<* j  V  Z 


(2.  10) 


4>i  fc  0 


bix  ,  VZ 


(2.  10.  a) 


and  must  satisfy. 


+  M  +  k2  4>  -  ° 


(2.  11) 


The  following  boundary  conditions  on  (J)  ^  and  hence  on  may  be  ob¬ 
tained  by  recalling  (2.  b;  to  (2.  9)  and  the  fact  that  the  walls  of  the  structure  are 
perfect  electric  conductors: 

al  ^t  =  0  it  a  =  o,vz  =>  il  =  o  at  x  -  o  v  z 
a)  ST  J  a* 


b)  =0  at  ^  =b  Z>o  ^  ii  =c  it  y 

6x  I* 


z>  o 


c)  continuous  a'  x  b,  t  <  0 


-ikz 


^(b-ro,  2)  -  </>  fb-0,z)  =  e  ;  Z  4. 


where 


(fltbi-o,  i)  •  <X,Z^*  = 


X  =  b-o 


X*  b+s  =  b-t  j 


where  tf  is  arbitrarily  small 


d)  ^  continuous  at  x  b  V  Z  "P  ■»  o,  2  )  rr  ^4*^  °j 

ax  J21  *-*• 


where 


(.  tj-t  o 


Ixs  bt  0 


e  and  hence  must  be  decaying  waves  since  a  loss  has  been  in- 


eluded  as  y  sCa^tZ2-)  ~ ►  *>  for  X  S  b 

Z  <  0  f  o£x^fc>* 


and  for 


and  1 7(^1  satisfies  the  edge  condition,  Meixner  [1954],  as  the 

edge  of  the  plate  is  approached,  that  is,  (j)  goes  to  zero  as 

\.'2  ,  -1/2  _  . 
z  and  |  7^  I  as  *  for  x  =  b,  z-0  -  J  ,  /-►  o  . 


Define 


x  =  -t-  z  T 


(2.  12) 


eo 


$(xa)  =  j  <F*;z)  e  *  dz 


(2.  13) 


*jk,  J  4>(x,z)  eJIZJz 


fir1 


'“y^J  ^'z)e 


1<XZ 


d  z. 


(2.  14) 


(2.  1  5) 


—  oc 


I  i  , 

£(*.*)  *  lcx.ee)  t  $Cx,<0  *  i-  i 

‘  +  ~  '  Tz^.i  lx 


1X2. 

e  d  z 


(2.  16) 


From  ‘hr  behavior  of  (x,  z*  for  any  given  x  as  z  -►ao  it  can  be  deduced  that 

(|U/X.')  *s  ana.’v’ic  for  T  >  and  ^>(xta.)  is  analytic  for  TV  kv  Hence 
for  anv  given  x,  is  analytic  in  the  strip  -k^X  ^ ^ 

The  boundary  conditions  (except  e),  in  view  of  definitions  (2.  12)  to 


2  lt>i  miv  now  be  written  in  terms  of  transforms  as 
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a.  1)  ,*)  -  ° 

b.  1)  §Vb+©,*;  =  <$U-o  ja)«o 

c.  i)  |ck+o,t)  -  5_(b-o,«ti  =-  if~r~{jzw 

d.  1)  J(l +b;  a)  a 

£.1)  S(b,«J  ~  —  *  *or  r>-l,i 

T  ' 

and  JiVb^'c)  ^  £'/%  *s  X--+°°  ^ 

(Abelian  theorem  given  in  Noble  [1958]  ). 

-I  X&Z 

Multiplying  the  wave  equation  (2,  11)  by  (£T)  S  and  integrat¬ 

ing  from  — oo  to  ao  with  respect  to  z  gives 


with 


d  1  _  y2,  cc)  *  <3 

d 


y=  -i(^-cc2)Vi 


(2.  17) 


(2.  18) 


The  branch  cuts  used  in  (2.  18)  and  shown  in  Fig.  2  have  been  chosen  so 
that  yhas  a  positive  real  part  when  -k,<  T -C  k?.  General  solutions  for 

La  L* 

$<*,  X)  satisfying  (2.  17)  .  ii.  a  form  convenient  for  applying  the  boundary 


conditions  are 


fl(a)  C»sk(/*) '  C<x^b 


-  5(3-)  e  -t*  DC<l)  ^ 


(2.  19) 


{2.  20) 


The  solution  for  $(x,  y)  is  obtained,  once  the  unknown  functions  A,  B,  C, 
D  are  found,  by  using  the  Fourier  invt  rsn  integral 


i 


T  = 

Im  (a) 

y~  BRANCH  CUT  FOR 

a- PLANE 

!  Y  r  =  ja2-  k2 

C  k 

INVERSION  CONTOUR-^ 

a  =  Re  (a) 


^  -k 
\ 


>  >  i  > 

«  .  *  -  fc*  1  w 

C  hon  e  ot  braiu  h  ails  (or  Y  i  a  -  k  )  and  the  unitnur 
used  m  me  Fourier  inversion  integral. 
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<o  +  \T 


f  -2*  2 

t*7A,  ,r  ;2-21> 

The  fact  that  (|(x,a  )  is  analytic  in  the  strip  requires  that  the  path  of  integra¬ 
tion  be  in  th  strip,  that  is,  in  (2.  21)  -k^  <  7*<  k^. 

The  boundary  condition  (e)  and  the  fact  that  V  has  a  positive  real  part 
f°r  k^  requires  D(ct  )  to  be  set  to  zero.  Likewise  boundary  condition 


(a,  1)  requires  C(ot  )  to  be  zero.  Now'  at  x  =  b  we  may  write 
(£fb-o)  4"  <£(b -o)  =  ft  Cc<h(Yh) 


(2.  22) 


-  Yb 


$(b+c)  ■+■  a^to) »  &  ^ 


(2.23) 


5>b-c)  + yr  S;^(vy) 


(2.  24) 


I  I  ~  v  b 

(£  (b-f-o)  +■  $.(b+e)  S“Y  B  C 


(2.25) 


Note  that  the  a  in  the  argument  of  the  functions  in  the  above  equations  has  not 
been  written  for  convenience  and  will  not  be  if  it  does  not  lead  to  confusion. 
However,  ^  ,  A,  B,  and  *Y  are  still  functions  of  3L  „  Boundary  conditions 

i  i  i 

(b.  1)  and  (d.  1)  show  that  |(b-to)  =  J(b-o)  and  is  now  defined  as  |j(b) 


That  is, 


$(b+o)  =  ijb-o)  *  5jb) 


(2.  26) 


and  also 


1  5 


(^(b+o)  =  ^b-o)  =  o  (2.  26.  a) 

Using  these  facts  in  (2.  24)  and  (2.  25)  yields 

i  ~Yb 

J(b)=Y  ft  Si'nMvb)  =  “YB€  27) 

The  two  unknown  functions  A  and  B  are  obtained  by  solving  for  the  unknown 
lb)  . 


Define 

Q  =  JCb  +  o)  -  ^  fb-o)  (2.  28) 

+  -f-  4* 

which  is  obviously  analytic  for  7”  >  -k7.  Subtracting  (2.  22)  and  (2.  23)  and 
using  (2.  27),  (2.  28),  and  the  boundary  condition  (c.  1)  yields 


where 


Z 


YztP  (a  -  b.) 


_ $j>) _ 

b  (a  +  k)(a.-k) 


(2.  29) 


L(<0=  — - SlhyL~ —  (2o30) 

The  function  L  has  branch  points  at  k  and  -k  and  {refer  to  the  definition 
of  Y)  L  is  analytic  for  -k^  <  k^  Therefore,  L  may  be  factored  into  the 
product  <L  L  which  holds  in  the  strip.  The  function  L^_  is  analytic  for  T >  -k^ 
and  L  is  analytic  for  7'<  k  .  This  factorization  is  the  difficult  step  in  this 
class  of  problems  and  it  is  obtained  by  a  limiting  procedure  as  discussed  in  sec¬ 
tion  2.  3. 


Let 
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with 


E(i)  =  E(«hLW  = 


Eja)  =  j  z  h 

iziP  (CL-k) 


j  fa+k)  Lfr> 

-(z/P  (a  -  fe) 


;  analytic  for  7"  < 


(2.  31) 


( 2.  32) 


£,(ol)  = 


^(a+fe)Lte)-2kL^)];  analytic  fo 


r  7">-k. 


(2.  33) 


1~z¥\  (a-k)  L'  — - - -  -  *2 

Multiplying  (2.  29)  by  (L+  (a  )  ).  (at  k)  and  using  (2.  311,  (2.  32},  and  (2.  33) 


gives 


D+  (a+kl  LW -EW  =  E>>  -  b  LtaT 


(2.  34) 


which  holds  in  the  strip  -k^  <  7"  <  k^.  The  left  side  of  <2.  M  is  analytic  for 

and  th«  right  side  is  analytic  for  T  k  .  Therefore,  one  side  is 

the  analytic  continuation  of  the  other  and  they  may  both  be  equa’ed  +o  a  function 
J  (  a)  which  is  analytic  over  the  whole  OL  *  plane. 


D+  ktk)  LW  -  E.£>  =  Ej«>  -  blA)W  =  J(a) 


(2,  35) 


The  application  of  the  edge  condition  gives  thp  unknown  function  J  (  C.  h  The  fac 

-1/2 

torization  obtained  in  section  2.  3  shows  that  L  la  j  behaves  as  Q. 


for 


- 1/2 


a-*  oo  ,  7*  >  -k^  and  L  ( £L )  behaves  as  a,  '  as  a-*cc  .  7~  <L  k^.  The  edge 


-  1 


i 


condition  (f.  1)  shows  that  D+  behaves  as  a  as  a.-*-  ao  .  7“  >■  -  k-,  and  ijb) 

behaves  as  a  ^  as  CL-+co ,  7~<  k?.  The  definitions  of  E  and  E  given  in  (2,  32i 

c  ~  ^ 
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-1 


and 


and  (2.33),  respectively,  show  that  E  behaves  as  cL  as  <£-+■  <x>  ,  7"<  k^ 
t  .  -1/2 

E  behaves  as  X  as  a.~*oo  ,  7“  >  -k  .  Therefore  J ( <X  ) ,  which  is  analy-- 

T  C* 

-  1  /2 

tic  in  the  whole  of  the  X  -plane,  behaves  as  a  as  a.-»oo  ,  7~<  k  and 
as  Ct  as  oo  ,  '7'  >  -k  .  Hence  from,  the  extended  form  of  Louiville's 

Cd 

theorem,  for  example  Hille  [1959],  J(a)  must  be  identically  zero,  that  is 

|'(b)  *  (a-fe)  L £a)  (2.36) 

Use  of  (2.  32)  and  (2.  27)  gives 


A  (JT«0  =  fl(a)  CosU^x)  *  2  LgO  Co5<i 


’I") 


"|  2- 77  Y  «Sihlo(rb) 


(2.37) 


B.(z,a)-BCaleYXS- 

TZ'fr’  Y 


Yb  —  Yac 


V  ) 


(2.38) 


The  formal  solution  is  obtained  by  inverting  (2.  19)  and  (2.  20)  and  may  be  writ¬ 
ten  as 


co+jr 


-ixz 


^LJ(x'a)e  °"x<-b 


(2.39) 


(j)(x;z) 


-  ~-lG~  z 

SJ_  \  B/sCjCL)  ^  dd.  )  ^£X 

fzr»-Ja>+lr 


(2.40) 


with  -k  <  7*  <  k  in  (2.  39)  and  (2.  40). 

w  w 
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2.  2  Choice  of  a  Closed-Region  Structure 

The  factorization  of  L(fl-),  given  by  (2.  30),  is  obtained  in  section  2.  3 
by  taking  the  limit,  as  a  parameter  approaches  infinity,  of  the  function  anti  fac¬ 
torization  which  occurs  in  a  related  C-R  problem.  The  C-R  structure  should 
yield  a  function  whose  factorization  is  obtainable  and  the  limit  of  the  function, 
as  the  transverse  dimension  approaches  infinity,  should  be  the  O-R  function 
(2.  30). 

It  should  be  pointed  out  that  it  is  not  necessary  to  resort  to  a  C-R 
structure  to  accomplish  the  factorization.  It  is  certainly  possible  to  intuitive¬ 
ly  choose  a  function  which  may  be  factored  and  be  such  that  some  sort  of  math¬ 
ematical  limiting  process  on  the  function  will  yield  the  O-R  function.  The 
choice  of  a  C-R  structure  seems  to  be  the  easier  choice  to  make  and  has  the 
added  advantage  that  all  mathematical  results  must  be  consistent  with  the  phys¬ 
ical  phenomenon. 

The  chosen  C-R  structure  for  this  problem  is  shown  in  Fig.  3.  The 
choice  of  a  related  C-R  structure  is  not  unique.  The  one  chosen  here  is  con¬ 
sidered  to  be  the  obvious  one,  that  is,  the  simplest  way  to  convert  the  O-R 
structure  to  a  C-R  structure.  The  O-R  structure  is  obtained  from  the  C-R 
structure  by  the  limit  of  a,  c  — ►  while  maintaining  a  -  c  -  b. 

The  solution  for  the  electromagnetic  fields  inside  the  C-R  structure 
is  formulated  in  the  same  manner  as  for  the  O-R  structure  in  section  2.  1.  The 
only  major  change  is  in  the  boundary  condition  (e)  which  must  now  be 
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re  corresponding  to  Fig.  1. 
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e)  6  $t  —  o  ht  x.  -  o.  .  V  Z  =>  ^  -  o  »£  2:  =  V  Z 
d  x  ; 

or  in  terms  of  transforms 

1 

e.  1)  $  (<x  J  a)  —  O 

This  change  will  be  reflected  in  the  general  solution  for  (x,  CL  )  in  x  ^  b, 
that  is,  (2.  20)  will  now  become 

{&(*;<*■)  =  B(a)  Cosh  D(cl)  5  » r?  h  \><X<0*.y2.4\) 

Duplicating  the  arguments  used  in  section  2,  1  will  'hen  give  the  following  re¬ 
sults: 

C(ci)  =  D^)  =  C  (2.  42  I 


k-.  m  sec  tion  2.  i,  The 
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-  1/2 


as  QL  oo 


asymptotic  behavior  for  K^(<z.)  is  chosen  such  that  K^(CL)^  CL 
for  7~  >  -k  and  K  (a  )  > — ■  <x  ^  as  cl— ►on  for  "7“  <.  k  .  With  these  condi- 

fa  “ 

tions  the  formal  solution  for  the  C-R  structure  is  given  by 


cr  ft  7 


-zaz 


(6(x,/)=  J_  ^  (L.  Ja  ;  o  ^  b 


yZT'  _ 


J 

OO  +  ZT 


(2.47) 


a>+t  /  ■  ^  ^ 

4>fe,z)  =  _t_  ^  D, (*,*•)  €  da. 

'  -CO  *ht  7" 


b  ^  a  ^  cx 


(2.48) 


with  -k2<7"<C  (2*4^)  and  (2.47). 


2.  3  Factorization  Obtained  by  a  Limiting  Procedure  on  the  Function  Appro¬ 
priate  to  the  Closed-Region  Problem 

A  formal  solution  for  the  EM- fields  of  the  structure  shown  in  Fig.  1 

is  obtained  once  L  and  L  of  (2.  30)  are  found.  This  factorization  will  be  ob- 

+ 

tained  by  talking  the  limit  of  K(CL),  given  in  (2.  45),  and  its  factorization.  Re¬ 
call  that  we  are  interested  in  a  factorization  using  functions  whose  common  re¬ 
gion  of  analvticity  is  the  strip  -k^  <  7  (Im  d.  )<  k^  and  whose  product  is  L(CL  ) 
for  values  of  d  within  the  strip.  This  strip  of  analyticitv  for  L(  CL  )  is  also  the 
strip  of  analyticity  for  K(CL  ).  Also  recall  that  Y  as  defined  in  (2.  18)  has  a 
positive  real  part  for  CL  within  the  strip.  Hence,  for  any  CL  of  interest 
(within  the  strip),  the  limit,  as  suggested  by  the  physics  of  the  problem,  gives 
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Lino  K(aj  =  L  i  no 

•X}C  — ►  CP  CL,  C  — *■  sc 

a.  -  c  =  b  ol-c  *b 


SjptYYfcL  Sinio(Yc) 

V  Si  n  io  (Ytfw) 


Li  i-r  K(a^  = 

CL  ;  C  ~>c to 
a  -  c  =  b 


-Yb 

e  =  b  Lw 


(2  49 j 


Li-*  Su  krL)  --  Li i-i 


a..c-»'Oj  c 

a.-e-b  ''>l 


L/ol) 


a? 

JL-  C  =  b 


"V 

L  (c^-fiTDc  -(<j'x-hi7i)c 

)  e _ =_e _ 

)  (<*z  -(.<fz+i'7Z)<x 

v  e  -  e 


— Yb 

=  e 


(2.  SO 


b>.Ca;s>.  Y*  &2-rt7^_  with  CJz  i>  O  for  any  <3L  within  the  strip.  Thus,  ’he 
‘uncuor.  to  If  far’  >r  d  for  the  O-R  problem  is  the  limit,  as  given  by  (2.491  of 
I’l  icut  ’i  to  b-'  bi-;*or«d  for  the  C-R  problem.  This  fact  enables  us  to  obtain 
•  a-:*oi  ;  1  .i,  &- ,  bv  *ui?  Lmihng  procedure. 

Th-  ■  nation  of  K(  0L  ‘ ,  given  m  (2,45),  may  be  realized  bv  using 

i‘*.  t-ac’  at  K(CL*  i"  a  ra’io  of  in’egra.1  functions.  Hence,  K(CL)  may  be  «.  x - 


p  r  «■  s  ^  »•  i  m  *  o  :  p 


-nhtu’p  products  from  which  K^(CL)  and  K  (CL)  may  be  t  on 


.pi,;  *.!'.  .•>.«>•  a,  ■  <i.  The  results  ar*» 


—  l 

K (<L)  =  Siw-^b)  S.n*bc)\ 
*  [k  Sin  (ko-)  ! 


_ab  r  T 

<o  cl  TT/ij-a\<^.7M  -^(d) 

ir,(i+s)e  jLir.d-w-Je  e 


_  iL 


TTd  +  ^-le  ** 

/»=»  <2-n 


Kjcn  =  K  t-aO 


(2.  51' 


2.  -2 


wi'.i  r  t  • : 
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(2.53) 


(2.  54) 


The  exponential  factors  are  inserted  in  the  infinite  products  in  order  to  make 

them  converge  uniformly.  The  function  ^  (CL)  is  chosen  to  be  analytic  and  to 

- 1  /2 

make  K  (CL)  behave  as  CL  as  qL -+cd  for  T"  >  -k  as  required  in  order  to 
4  2 

be  able  to  solve  the  Wiener-Hopf  equation  for  the  C-R  problem.  The  function 
(cL)  is  obtained  from  a  knowledge  of  the  asymptotic  form  of  K^(CL).  The 
asymptotic  form  of  the  infinite  products  is  obtained  by  comparing  them  with 
the  infinite  product 


«  . 

TT(i  + 


q 

an  +  b 


g. 

an 


rtte+0 

p(* 


*2.  54.  A) 


where  P  (a.)  is  the  Gamma  function  and  p  is  Euler's  constant.  The  asymp¬ 
totic  form  of  (2,  54.  a)  is  obtained  by  the  use  of  Sterling's  formula.  Since  each 
infinite  product  in  equation  (2.  52)  behaves  in  the  manner  of 

=  i  ftTT  -h  CKn”')  (2.  54.  fc 

for  large  n,  their  asymptotic  form  may  be  obtained  from  (2.  54.  a)  to  yield 


tic L)  =  _  i-i  1C  -t  b  in  l°t\ 


t  (-  +4 - 


imr 


rit 


<2.  54.  c ) 


The  function  X  ( CL  )  results  in  K  la)  behaving  algebraically  for  the 
'  + 


solution  of  the  C-R  problem.  However,  when  merely  using  K(q)  to  obtain  the 


0-R  factorization  it  is  not  necessary  to  include  this  term.  Its  inclusion  does 
lead  to  the  possibility  of  obtaining  the  O-R  solution  from  the  C-R  solution  as 
'.v ill  be  discussed  in  section  2.  5. 


K  t  ( £t )  and  hence  K  ( <X  ),  without  in  luding  the  X  (  &  )  term,  may  be 


expressed  in  an  alternate  form.  This  will  prove  to  be  advantageous  when  lim¬ 
its  are  to  b*»  taken.  In  particular,  ont  can  write 


K_C&)  *5 


*  , 

Si4b)  •  7T  ( 1 +  ) e* 

k  -^OecO,  !  n-  ^ 


l 

Ht<0 


(2.  55' 


with 


where: 


HCa)  —  Li*»  ^  aA  did 

-'7ri  i  ; 


“FfCL.  u) 


(Kw1)  *  Fi(^  _ 

Fit")  ^u») 


(2.  5bi 


(2. 


(2.  58 


a)  F  (-)  has  simple  zeros  at  JoJT  *,  H  *  I  .  2.  ,  •  •  • 

1  c  '  J  * 

b)  F  (•-’)  has  simple  zeros  at  Q3T  ]  r  ~  1 )  2- ,  *  *  • 

z  <x  ' 

c)  F  U)  =  .  F  '(v)  - 

1  doo  2  Ju, 

d)  ^  in  the  contour  shown  <n  Fig.  4  with  o  <Y  TT/o^  £  <.  ^ 


Ftg.  4  Contour  used  in  the  representation  of  K  (  a  )  given  by  bb). 
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NT  <  <j  <  (  N  t  £  T  Is) 


The  functions  F  (to)  and  F  (to),  may  be  identified  with  the  functions 

1  w 

which,  when  set  to  zero,  result  in  the  characteristic  equations  which  give  the 
transverse  wave  members  for  the  regions  b<L  x  ^  a,  z>o  and  o  <  x  <  a, 
z  <  o.  respectively,  in  the  C-R  structure,  Fig.  3.  In  this  particular  prob- 


F^  (to)  =  Sin  (ci.c) 
F^(oo)  -  Sin  (toa) 


(2.  59) 
(2.  60) 


That  K,  (  )i  given  by  (2.  55),  is  the  same  as  K  {CL),  given  by  (2.  51) 
may  be  verified  by  integrating  (2.  56).  For  any  QL  with  7~>  -k^  the  integrand 
in  (2,  56"  is  analytic  with  respect  to  to  in  -k^  Im  to  <  k^  except  for  the  simple 

i  i 

poles  at  the  zeros  of  F  (to)  and  F  (to).  In  this  case  F  (to)  and  F  (to)  do  not 

1  U  Yu 

have  any  poles  within  ^  ,  Using  the  calculus  of  residues  one  obtains  equal¬ 
ity  between  (2.  51)  and  (2.  55). 

The  expression  for  H  (  &  )  may  be  expanded  further  by  integrating 
along  the  path  of  integration.  That  is 


d+ze  0+t€ 


,  -  4  6 


.  n  n. 

d-tfe  d+16  orl-S 


-z  ( ■  -y) 


(2.61) 


Now 


'■■*Y Nbj  means  the  largest  integer  in  Nb  . 
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o  +  i($-r)  o- z€ 

s  -t  s 

oti  e  o-ice-r) 

because  the  integrand  is  an  odd  function  of  to. 


Also 


4 =  -vi 


(2.62) 


(:.  63) 


d  +  ife 

d-i6  N 

,  -2 

because  the  integrand  behaves  as  to  for  large  to. 

Therefore 


(2.  64) 


<33 


(2.  65) 


where  £  is  any  positive  number  such  that  0  <  6  <1  k  ,  that  is,  H  (d.)  is  inde- 

6 

pendent  of  £  .  A  useful  representation  for  K(GL )  is  now  available  and  may  be 


written  as 


KU)  -  ■Sin In  (Vb^  S ih(kc\_ 

Y  Si  h( ko.) 


C  Hi(a')  •+-  H  (-a)"] 


(2.  66) 


with  H  (cC  )  given  by  (2.  65). 


A  representation  for  L(C£.)  for  any  <x  within  the  strip,  -k^  <  T k^, 
ana  from  which  L+  and  L  may  be  chosen,  is  obtained  by  taking  the  limit  of 
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\Z.  66)  as  a,  c  -►  co  while  maintaining  a-c-b.  The  limit  of  the  left  side  oi 
(2.66)  for  values  of  ex  within  the  strip  is  given  by  (2.49),  that  is,  bL( CL  ), 
Now 

i.k  b 

Lim  s.r  (Itc)  =  e 

9--  C  =  b 

and  let 

H  (tf.)  =  Li* 

(X  c-»oo 

ol- t  •  b 

The  function  f(  <X  ,  co'*  is  independent  of  a,  c  and  the  function 
this  problem,  is 


12.  67; 


•  h  68 


g(J?i  *  'h  for 


oC|  +  Z€)=  C  CosfcC*- Ticji 

*  SlhCcC5+l€)] 


a  Cos 

5|n  [jxCf  + 


[Z.  69) 


Hence 


l ii m 

0-,C  cw 
a  -  c  -  b 
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( 1  :p  i£  )  = 


(2.  70, 


Define 


&(l) 


L‘im  q(%-l€.)  ~  L\h-  q(j  +  i€)l! 
S.'tTZ  j  a.(c-*-oo  J  Q-jC— ►  oo  J 


(__a.-c  =  b 


cl—  c  =  b 


( 2.  70.  a) 


I  £-*■© 


then 


OCP 

HCcl)  -  J  ft  <*-,&>)  d  ^ 


2*  71) 


60  =  O 


1 


2  9 


where 


<  2 ,  7 1 .  a  1 


It  is  interesting  to  observe  the  physical  significance  of  it.  ”0.  ah 
G*(to)  given  by  (2.  70.  a)  is  equal  to  the  difference  of  A^  and  A  ,.  A^  a^d  A., 
are  the  inverses  of  the  spacing  of  fhe  zeros  of  Fj(co'i  and  „  Tvis  is  ex¬ 

pected  if  one  had  worked  with  the  inf  .nde  product  r'orm  o'  K(CLh  The  natural 
logarithm  of  K(  CO  would  convert  products  over  n  to  sums  of  logarithms  o"-‘r  n, 


One  could  then  change  the  summation  to  a  sum  o‘-er  ntr /  a  and  ntr/c,  that  is.  *he 

zeros  of  F  (co.  and  F  (co),  As  a  and  c  approach  infinity,  iur  =  and  ou*  = 

1  ,  I 

would  approach  continuous  variables  w  and  co  „  Also,  Au  v/a  and  Ao 

n  n 

tr/c  would  approach  <io  and  cko  „  In  the  lirn:4  one  would  ob’a’n 

oo 


<L)U J) 


OJ-O 


-L-I\  <U 

A  O0y\ 


1  -7  ) 

C-  -  • 


but 


b 


I  am  IT  ( 2.  73 1 

Comparing  (2.  7 2 1  with  (2.  71  ■  justifies  the  interpretation  of  G  ko)» 

This  method  could  have  be«n  used  but  if  is  not  as  convenient  or  as  suf¬ 
ficiently  general  as  the  repr  es-:ntaTion  given  bv  equation  \2.  5 h  when  analyzing 
more  difficult  problems  as  in  sections  ^  and  4,  However  it  does  provide  a 
check  on  the  limit  of  g(£-  h  The  dif^reno  of  the  in’  ”rsc  of  spacing  of 
the  zeros  of  F  (co)  and  P'^fcoi  can  be  calculated  for  increasing  ’  alues  of  a,  c 
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a-c=b,  and  the  result  should  approach  the  expression  obtained  for  G  (oi), 

(2.  70.  a).  The  check  in  this  case  is  trivial  as  the  zeros  of  F^(u))  and  are 

obvious. 

The  limit  of  (2.  66)  and  the  results  given  in  (2.  49),  (2.  67)s  and  (2.71) 
show  that  L(0L )  is 


LCa)  = 


Sin  1i(yL>) 

Yb 


ikL  Lh(«.hh(-cl)] 
e  *  e 


(2.  74) 


with  H(O-)  given  by  (2,  71).  The  choice  of  L^(CL)  may  now  be  made  as 


ii!>  +H(<0 
e  * 


(2.  75) 


*  L(r<*-)  (2.76) 

f 

The  range  of  analvticity  of  L  (QL)  can  be  deduced  from  the  regions  of  analytic! - 
ty  of  the  infinite  product  and  H[CL)»  The  imaginary  part  of  3n,  (2.  51),  is  great¬ 
er  than  k  and  hence  the  infinite  product  in  (2.  75)  is  analytic  for  '7'  >  -k  . 

L  “ 

The  function  f(  CL  ,  co )»  (2.  *7),  satisfies  the  conditions  of  theorem  A,  page  1  1, 
Noble  [1958]  and  also  restated  in  the  Appendix.  Hence  H(cL)  is  analytic  for 
7^^-k^.  Therefore.,  L,(CL!  is  analytic  for  *7">  -k^  and  L  (OL)  is  analytic 
for  7"  ^2*  t^ie  desired  range  of  analyticity.  The  product  of  L^(OL  )  and 

L  (CL)  is  obviously  L(CL). 

The  functions  L^(CL)  and  L  (cl)  can  be  arranged  in  a  slightly  differ  nt 
form  which  is  preferable  for  numerical  calculations. 


•'2,  r ..  r  vL1  PriGcD-  HC-aD 

LW-  [Lw  L«a  .  jf^l+^e  .  SOW)  eY  .  e 

+  L  Ca .)]  „  a.  ig  Y  ? 

L  '  J  -rr('-t)e^  L  J 

ns  -J 


.  ?  ^7  . 


!r  •>., *  section  equation  *2.  77;  will  not  be  us«rl  but  >ts  'om  *v  ]]  be  used  m  see 


'ion  1, 


»: at  .or.  [I.  71)  may  be  r  egr ated  to  obtain  H(  CL'-  and  hence  I.  ,ol*  ;n 
u  Th.a  n-i'- jr«  ot  H(flL)  -,n  •qna'.roti  (2.  71;  m-»v  v  e  o'?'a  one.  t*  • 


lof  ‘ d  form. 


.fo’ lo-wne  i:> ‘ ea  r  ^ t; on  :s  p^tormed. 


XE"(i+c^r)'/i) 


* _  JcO 


Con  si  <let  Cun  cm  or  O  di  w! 


*,av“>  =  to  lB(|  +  _  2  _ 


Z'!'z  I  f  .  ^-a*n9lf-<fl-w) 
(K^_a;  + a. 


T>>e  branch  o.t  In  (k  jj)  :s  chosen  to  be  in  the  lower  halr  plar  *>  and 
branch  of  1^  <k  oj'  if  chosen  to  S~  m  the  upper  half  plane.  lf^ncr  for  anv  CL 
suer:  rha  •  7“>-k?J  Q  CL .  w'.  is  an  ana.’vtic  (unction  o'  u>  'or  ‘^e  lm 
r u t;  ,1,-nvaU’-?  of  Q  a  ,  w.  wj»h  respect  *ou>  for  anv  w  :n  *b;s  rant»**  -s 
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Therefore,  the  integral  (2.  77.  a)  becomes 


00 


Q(^w) 


--U+ilfe-lY  !„(*-*.) 
2  V  *  / 


(2.  77.  c, 


and  hence 


H(<L)-  iby  I«(sl£3C\  +itfe  -  ihk 

nr  \  k  /  ^  2 


(2.  ?8) 


and 


•*.  -  ~  tab  .  ib_  V  L/.g-OC 

L(a) = Is^W® , + i)  e1^] .  e  ^  ^ 


-1-XCa.) 


(2.  79) 


With  L_(a)  =  L+(-a)  and  X  (- «-)  =  Oc  (GL).  The  principal  determination  of 
the  logarithm  is  used,  that  is.  In  (1)  =  0.  The  unknown  factor,  exp.  \^-)C  (CL-TJ s 
which  must  be  analytic  in  '7*>  -k^,  is  added  so  that  the  asymptotic  form  of 
L+'a)  is  algebraic.  The  algebraic  behavior  of  L+(C-),  and  hence  L  (CL),  en¬ 
sures  algebraic  behavior  of  J(CL)  and  hence  its  determination  in  equation  (2.  3 5 • 
by  the  extended  form  of  Liouville's  theorem. 

For  large  CL  in  7">  -k^,  we  have 


(2.  80 ' 


The  asymptotic  form  of  the  infinite  product  may  be  found  by  the  use  of  equa¬ 
tion  (2.  54.  a)  and  gives 


TRh-JH  e 

r  « t 


ab 

tnTT 


CL 


(2.  81' 


X 


3  3 


where  P  is  Euler's  constant.  Therefore  if  X  (CL)  is  chosen  as 


X(<o  * 


then 

Lta)  - 


‘/2 

Sin(leb)  * 
feb 


(2.  82 


,  l  /? 

with  L1*.  a.  ^  a.  as  GL-^CO  for  7"  >  -k^.  The  function  L  (OL* 

,  -1/2 

Wi*h  L  'a  ^  OL  '  as  a.-*00  for  7*<.  k^. 


f  <\  8  3 


CL 


2.  4  Sol’. Mon  of  *he  Problem 

Tve  factorization  given  by  equation  (2.  83)  and  the  transformed  field 
quantities  g:'-en  by  equations  (2.  37),  (2.  38  ,  (2,  39),  and  (2.40'  permit  the  de- 
•erm'natmn  of  the  field  quantities  of  interest.  The  scattered  fields  wj*hin  the 
waveguide  are  given  by 


oo -e  t  r 


L Ja)  Cosh  C'fx) 

y  Sf*k  (Yt) 


-i<LZ 

0  dcL  * 


Z  >o 


(2.  84) 


wi»h  -k  >  <  7-  <  k,. 

V.  tj- 


Since  i  >  0  and  the  integrand  has  no  branch  cuts  in  the 
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lower  half  plane,  the  path  of  integration  may  be  closed  in  the  lower  half  plane. 
The  solution  for  ♦  is  then  obtained  by  using  the  calculus  of  residues,  for  ex¬ 
ample,  the  reflection  coefficient,  R,  in  the  waveguide  is  given  by  the  residue 
at  0[=  -k. 


R  =  -li  <k)  =  -l/  (k  ) 

+  TO 


(2.85) 


Higher  order  reflected  modes  are  available  by  the  same  procedure.  Note 

that  once  a  result  as  given  in  equation  (2.  85),  or  those  following  in  this  section, 

are  reached  k  may  be  taken  as  real  and  equal  to  k^.  This  may  be  done  because 

k  is  infitesimally  small  compared  to  k  and  the  evaluation  of  a  function  at  k 
Z  1 

is  equivalent  to  its  evaluation  at  -  co  =  *• 

The  fields  for  x>  b  and  hence  the  radiated  field  is  given  by  equation 
(2.40)  with  Bj  (  x.j  CL  )  given  by  equation  (2.  38).  The  function  B^x,  c*.  )  has  a 
branch  cut  in  both  half  planes  and  he  .ce  little  is  gained  by  closing  the  contour 
by  an  infinite  semicircle.  The  radiated  fields,  asymptotic  behavior  of  <p  given 
by  equation  (2.40),  is  usually  obtained  by  an  integration  procedure  known  as 
the  saddle-point  method,  for  example,  Morse  and  Feshback  [1953].  However, 
in  this  class  of  problems  the  use  of  equivalence  theorems  for  fields,  for  ex¬ 
ample,  Deschamps  [1962],  will  prove  to  be  more  direct  and  convenient.  The 
field  in  the  aperture  is  given  by 


r  -iaz 

\  B,(b.a.)  e  da.  -  VZ 


i  > 


fiib+O/Z) 

Recall  that  H  (b  +  o,  z)  -  6  (b  t  o,  z)  -  d)(b  +  o,  z). 

v  '  t 


— cp  +  \T 


(2.  86) 

Therefore,  H  m 

y 
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the  aperture  is  given  by  equation  (2.86).  The  equivalent  electric  surface  cur¬ 


rent  in  the  aperture  is  given  by 

A  A 


*  * 


<X  —  ft  k  H  *  .2:  x,  ^  *  Z  (y+Oj  Z)  Six— y)  (2. 87) 

The  actual  far  field  is  equal  to  the  far  field  due  to  the  equivalent  electric  cur¬ 
rent  2J.  The  vector  potential  due  to  the  current  2J  is  obtained  by  knowing 
the  asymptotic  form  of  the  Hankel  function  ■  as  described  in  Har¬ 

rington  [1961]  and  gives 


iko/* 


CD  00 

n  n 
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z 


1  ~lZ7T 


- 1  a  •  ^ 

(jKb+o^  Zj)  6(x1-y)  e  Jx.L 


7°  z1--^  -L 


CP  X.--CP 


(2.  88 


where: 


k  «k.£Vste)z  +kc  s‘tn(e)x 


rx  =  zA  z+j:^ 


I  A 


/ 


=  (x^-hZ2)1'1 


Integrating  (2.  88)  with  respect  to  x  gives 


/?=  e 


!  “  -jbjZ 

il^Y 


«ttT 


2  7T 


-1  T  -ZCLZj.  x 

^  J  B(Cb;a)C  d<xj  dZj 


5-00 


and  hence 


(2.  89, 
1.2.  90 
(2.  91' 


!.  92] 


c$  (x'  is  ti  e  delta  function. 
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-r  a 

R  =  t- 


-fly 


—  Z  -Sih  (d) 
1 - 


B>,-kc0s©) 


From  Maxwell1®  equa’ions  w«  ob*a>n 


(2.  Q3) 


H  =  -ik.  FI  Sir  (6) 


or 


Sin (6) 

» 


(  2.  94 


In  this  problem  H  in  the  far  ft. eld  is 


3  l'lff* 

with  0  <1  0  *C  rr. 


CoS  s) 


(2.  9  5; 


The  real  power  reflected  in  *he  w’aveguide  and  »he  real  power  radiated 
in  the  space  wave  arc  obtainable  once  »he  number  of  modes  propagating  in  the 
waveguide  are  known.  If  th«  guide  dimension  and  fn*q’i«‘ncv  of  ope*afion  aie 
chosen  so  thaf  only  »he  TE.M  mode  can  propagate  the  resuPs  are; 

Normalized  Keflec'ed  Pow.-r  _  |r|  c  1 1.  »k  |  if!.  Po  i 

The  normalised  power  radn*.  d  :s  go  en  by  Povn'mu  '*•  c'or. 


Normalized  Rad.a'»d 


sin  e  the  inc  ieN  n*  po'* 


Power  W 

r  a  d 


•risk  b  _  £ 
o  < 


o  _  _ ♦  ^ 

a-  \  £x  H  •  o  p  d  8 

bob  J  I  ' 

e*~ 


\  cw 


i  f  .  <  ”  i 


*w 

§■0 


Wrs<)  ■=  X 


In  this  problem  (2.  98)  becomes 

t r 


Ccs  $)  B 


Ae 


W^-AfejL^wf  j  lL(h.C..  «)f  de 

The  radiation  pattern  as  a  function  of  0  ic  giver,  by 

I  Hjl  «  |B,U,-I».C0j  »)  •  St*  «| 

i  nu  in  this  problem  is  proportional  to 


2.  98) 


(2.  99) 


\i.  1-  ) 


Hjjl*  1  Cos  b)\ 


(2.  101) 


The  numerical  value  for  W  ,  ,  W  ,  ind  the  far  fic-i  oatterns  are 

ref.  rau. 

given  in  section  4.  Some  of  the  equaticns  in  this  section  have  been  obtained, 
with  sufficient  genet ality  so  that  they  will  apply  in  the  analytic  of  the  problems 
in  section  3. 

In  this  section  1,^(0.)  and  hence  L  (ct)  were  available  in  closed  form. 

Hovx/ir,  this  is  not  always  feasible  as  the  integration  usually  cannot  be  carried 

.  t_  1  ce  \  ( OL )  i*  obtained  in  closed  form,  for  example  (2.  83),  the  os*  may 
+ 


>c  'i-rauced  to  zero,  betting  <3  to  aero  results  in  k^  =  0,  k^  =  i 


n  -  w 
o 


?o  *0 


3 


g'  =  o,  k  -  k  in  (2.  83).  The  path  of  the  Fourier  inver- 
v  1  o 


iu cfciore,  when 


38 


sion  integral  used  in  (2.  39}  and  (?..  40;  w.cen  is  greater  than  zero  is  shown 

in  Fig.  2.  This  path  must  he  indented  when  the  Joss  approaches  zero,  that  is, 

when  k — V  the  contour  is  shown  in  Fig.  5. 
o 

The  equation  for  L,,  ( a '  must  be  interpreted  in  a  simitar  manner  as 
the  Fourier  inversion  integral  wren  the  Joss  is  reduced  to  zero.  For  example, 

j  i 

H(  CL  )  used  in  (2.  75)  and  defined  bv  iZ.  7]>  becomes,  when  k- — ►  k  , 

o 


Htao 


(2.  102) 


with  the  contour  ^  given  by  Fig.  6.  However,  no  use  is  made  of  the  inte¬ 
gral  form  of  H(a.)  given  by  {2.  7 1 '  and  (2.  102)  as  the  closed  form  is  available. 

In  the  more  difficult  problems,  discussed  m  trie  following  sections,  the  integral 
form  is  used  and  the  contours  of  integration  must  be  interpreted  as  discussed 


here. 


2.  5  Comments  on  the  Method 

The  factorization  obtained  ip.  seer  on  ?..  ?  /or  h(c£_)  cau  be  verified  by 
reference  to  the  problem  of  the  paraIJ<-l  plat.-.-  duct  discussed  by  Noble  [1958]. 
The  function  to  be  factored  ;s  the  same  n  eich.  problem.  A  comparison  of  the 
factorization  for  L,  given  bv  equal. on  ;7.  8'-),  with  Noble's  Solution,  shows 
that  thqy  are  the  same.  He  ohla.tr  -d  the  fa-'torizaHo  :  not  on  L(CL)  directly  but 
on  a  function  related  to  L(  cl  •  ar  d  t  •  use  c!  a  form-ii  factorization  procedure. 

The  form  o’  .L  obtained  by  means  of  a  Jim  Ting  procedure  on  K, 


«  <  . 
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(2.  66),  proved  advantageous  because  of  the  generality  of  H  (&.)  given  by  equa¬ 
tion  (2.  56).  In  the  problems  discussed  in  section  3  it  will  be  seen  that  g(w) 
changes  for  the  various  problems  and  not  f(  a. ,  o).  As  already  mentioned,  a 
check  on  the  limit  of  glj?  +  i£  )  is  available  because  it  is  equal  to  the  differ¬ 
ence  of  the  inverse  of  the  spacing  of  the  zeros  of  the  characteristic  equations 
F^(u)  and  F^(co),  The  form  of  h\o.),  equation  (2,  56),  also  applies  to  certain 
other  geometries  involving  coupled  waveguides.  It  is  found  that  in  these  cases 
the  form  of  (2.  56)  remains  unchanged.  One  needs  only  to  substitute  the  char¬ 
acteristic  expressions  for  F^(cr)  and  F^u)  appropriate  for  the  particular  geom¬ 
etry  under  consideration.  For  example,  this  method  of  solution  should  be  applica¬ 
ble  for  analyzing  radiation  from  a  semi-infinite  circular  waveguide.  The  above 
discussion  brings  out  a  strong  similarity  between  cylindrical  and  parallel  plate 
waveguide  problems  that  may  not  be  apparent  without  a  close  lock  at  these  prob¬ 
lems. 

The  form  of  given  by  equation  (2.  77)  will  also  prove  convenient 
when  numerical  calculations  are  attempted  in  section  4. 

The  C-R  problem  has  been  used  primarily  to  generate  the  O-R  factor¬ 
ization.  However,  the  relationship  between  these  two  problems  may  even  be 
made  more  general.  If  the  asymptotic  form  of  K^(flL),  namely  ^(<JL),  is  included 
in  (2.  55)  and  hence  in  the  limiting  process,  we  obtain 

Li*  K(a)  «-/?  LU) 

a  c—+oo  '+ 


(2.  103) 
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This  follows  from  the  fact  that  the  limit  of  the  asymptotic  form  o:  K^(  00  be- 
com<  -  the  asymptotic  form  of  L^(d.).  That  is,  for  K+(GL),  the  )C  (CL)  is 
given  by  (2.  54.  c),  namely 


The  series  involving  the  terms  and  ^  may  be  expressed  as 


(2.  104) 


f'M 

FM 


\OJ 


(2.  105) 


where  F  and  F  are  given  b^  equation  (2.  59)  and  (2.  60),  respectively.  The 

X  L 

contour  2  is  the  same  as  the  contour  X  shown  in  Fig.  4,  except  that  now 

w 

bd3T  <T  d  <  (m+|).ZL  .  Likewise  ^  is  given  by  Fig.  4  with  fc^UET  <  d 

<  (M  +  I)  TE  .  Integrating  along  the  contours  gives  (  &  is  set  equal  to  r) 


- - - 


H- 


+ 


l 

Z7T2* 


Mir- £  6  i 


v  • 

a )  =  o- 16 


[*rr  -(-  £  e 

"  J 

UziO-rie 


MjT-ie 
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Fj_ 
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cU 

Itf'-co'1-'  F, 

• 

'■f  7  t 

J 

(2.  106) 


Now  a  and  c — >■  ce>  but  so  does  M  such  that  M/a  — °''J  and  M/c  — 1 .  Hence 


we  may  write  as  a  and  c  become  large 

m  ,  .  r 


2.  (_L 


I  \  ^ 


it- 


M'jr/i 

r 


c\ 


u 

0}-  a 


Tr-fi?-  ^ • 1 


d  — 


C 


W  =  o 


_  A 


-\  Oi 


Hence 


2  ( -L  -  -1  \  =  -J-  [cl  In  ( ZMir)  -  c  In  ^  £  k  In  (-i  k) 

n--iUh  l(,  1  IT  L  V  cc'  )  \  C  /J  7T 


Therefore,  (2.  104)  becomes  as  a,  c 


( 2,  107) 


(2.  108/ 


LiVn  X  ^CC')  -  —  j.  OL  —  b  In  /-  Idkh)  -f  h  J  Liro  T  !h  (M  )  —  ^  '7T  i 

o^c.  — *  oo  7T  V  '~'7r/  r'-' 

a-  c  =  b  L  L  J 


( 2,  109) 

which  is  X  (fiL)  f°r  the  O-R  structure  as  can  be  seen  by  comparing  <2.  109) 
with  (2.  82)0 

The  fields  for  0  <  x  b  involves  the  Fourier  inversion  of  A^tx  QL  ) 
given  by  equation  <  2.  37)  for  the  O-R  structure,  and  of  C^(x  CL  ).  given  by  equa 
tion  (2,43)  for  the  C-R  structure.  For  any  gl  within  the  strip  and  any  0<x<b 
the  limit  as  a,  c  -*-ac  ,  while  a-c-b,  of  C  ^  (x,  cl  )  is  A  ^  (x,  jl  )  provided  (2,  1 0  3  i 
is  true,  Tha+  is,  K  +  (CL*  includes  ♦  he  term  X  (  JL  >  which  makes  ;t  algebraic 


in  the  proper  half-plane.  Hence  the  fields  in  0  <  x  <  b  for  the  C-R  structure 
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approach,  in  the  limit,  the  actual  fields  of  the  O-R  structure. 

This  establishes,  for  example,  that  the  reflection  coefficient  for  this 
chosen  C-R  structure  converges  to  the  value  of  the  reflection  coefficient  for 
the  O-R  structure.  Indeed,  this  is  what  Mittra  and  VanBlaricum  [1965]  re¬ 
ported. 

Likewise,  the  limit  of  D^(x,  a.  ).  given  in  (2.44),  for  b  <  x<  a  and 
Q-  within  the  strip  is  equal  to  B^(x,  ft. ),  given  in  (2.  38).  Therefore,  the 
fields  for  b<x<a  in  the  chosen  C-R  structure  approach,  in  the  limit,  the 
fields  for  the  O-R  structure  for  b  <  x  (a  goes  to  infinity  on  the  limit).  This 
establishes  th?t  the  O-R  solution  is  a  limit  point  of  the  C-R  solution  in  all 
space.  Therefore^O-R  field  quantities  may  be  obtained  from  the  correspond¬ 
ing  C-R  field  quantities  by  a  limiting  process  as  suggested  by  Talanov  [  1 9 S9  1 
and  Mittra,et.al.ri966]. 
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3  EXCITATION  OF  A  DIELECTRIC  SLAB  BY  MEANS 
OF  A  TRUNCATED  PARALLEL  PLATE  WAVEGUIDE 

3.  1  TE  Excitation  of  a  Surface  Wave  Structure 

The  excitation  of  a  dielectric  slab  by  means  of  a  parallel  plate  wave¬ 
guide  with  or.e  plate  truncated  is  analyzed  in  this  section.  T- e  surface  wave- 
structure  is  shown  in  Fig.  7.  The  relative  dielectric  constant  of  the  slab  <  jC  ) 
is  assumed  to  be  greater  fhan  one  and  hence  the  possibility  of  the  structure 
supporting  surface  waves.  The  incident  fi°ld  in  the  wav^gu.de  section  is  taken 
to  be  the  lowest  order  TE  mode  with  the  electric  field  intensity  parallel  To  the 
walls  of  the  guide.  A  slight  loss  due  to  finite  conductivities  (  (j  y  d  ?'  is  as¬ 
sumed  in  each  region.  However,  when  the  final  solution  ;s  obtained  *h:s  loss 
is  permitted  to  approach  zero. 

The  TE  polarization  is  described  in  detail  m  sections  3.  1.1  to  i  1.  4 
as  opposed  to  the  TEM  excitation.  The  details  of  the  formulation  of  the  prob¬ 
lem  for  the  TEM  excitation  follow  closely  those  in  section  Z  and  only  the  perti¬ 
nent  differences  and  results  are  given  ir  section  3.  Z. 

3.  1.  1  Fo.mulat-on  of  the  Problem 

The  tc*al  el  ctromagnetic  fields  are  obtained  by  solving  the  scalar 
wave  equa’ions  for  the  scattered  scalar  poTen:ial  .  De'me 

+  4  <5- i: 


4  7 


-  xS  ~~ 

C^  =  Sin^L*]  w  '  ) 


/ 


(3.2) 


(3.  U 


.  l/'2- 
\-  L  on ^  07  )  = 

Tc 


k,  f  <:  k  f 


(3. 


with  >0,  >  0.  Note  that  the  parameters  are  chosen  so  that  (2r/b)  > 

■>  ? 

Rek^  >(rr/b)  .  This  ensures  that  the  lowest  order  mode  propagates  in  the 
waveguide  and  that  the  dielectric  sU"'  is  excited.  The  scattered  scalar  poten¬ 
tial  tp  is  given  by  the  solution  of 


+  £A  +-  ^  t  =  ■ 

>  z3- 


il 

a  ^ 


b  ^  *• 


(  V  o ) 


The  constant  k  is  given  by  equation  (2.4).  The  electromagnetic  fields  are  ob¬ 
tained  from 


—  T 

-u  -f-t 


>>  t 


H  —  “ 


::  -ot>. 


i  i .  « > 
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H 


^  'K 


Z  itu/Jc  ^  (3.9- 

Defining  Fourier  transforms  of  cj>  by  equations  (2.  12)  to  (2.  16)  re¬ 
duces  the  problem  to  the  solution  of  ,  . 


cf  __  V2  <£(x;ol)  *  C  ) 

‘cl  ;c" 

J:'  1'.^^  __  yz  -  o  > 

•'  " 

where: 

2  ,  '/?„ 

y;=  (^-k, ) 

and  V  is  given  by  equation  (2.  18).  The  asymptotic  behavior  of  (x.3  zr  re¬ 
sults  in  (3)  (x,  cl  )  being  analytic,  for  7~>-M:n(k  ,  k  )  and  ^-(x.,cC-i 

3-i  2  4 

b^ing  analytic  for  7~  <  Min(k  k  ' 

c ,  4 

Reference  to  Fig.  7,  equations  (3.  7)  to  (3.  9),  and  equations  (2.  12)  to 
(2.  16)  gives  as  the  boundary  conditions  on  (£  (x;  cl) 

a)  {  ,  CL)  =  C 

b)  cL)~  'f/b -c;  cl)  =  ° 

c)  t  b +o; -O  —  (p_(b-OjCL)  ■=  C[  tb) 


0  <  JC  4  o 


b  4  oc 


(3.  ICi 


(3.  11) 


'3.  12: 


d) 


i_U  1  0f  CL)  - 


( b-e;  at) 


b 


A 


* Min  Minimum  value  of 


■x  -  <5/2  . 

e)  y)(b;cL)~a.  Sts  cl-*- oo  -for  T  <  Min  ( k* ,  ^4.) 

$(b,<0~  a  Z  as  CL—+"  oo  -for  |  7*  >  “"Min  (/ev  /*--*) 

Boundary  condition  (e)  reflects  the  fact  that  E  =»  (t)  and  X  ~  z+  a 

y  T  y 


at  x  =  b, 


A  solution  of  (3.  10)  and  (3.  11)  in  a  form  suitable  for  the  application 


of  the  boundary  conditions  is 


(jjU,a)=  /7(a)  Sir»fi(Y*)  +  C(a-)  CosV  (Y%)  )  o  <  ;*  _<  io 

_Yx  v  vC 

$(*,<0  -  Bftr)  e  ■+  DCo-'i  e  ;  * 


(3.  13) 


(3.  14) 


Recalling  that  Y  has  a  positive  real  part  for  any  a.  for  ^  and  the 

fact  that  we  are  seeking  decaying  waves  at  infinity  requires  that  D(  CL )  be  zero. 
Also  bo-undary  condition  (a)  requires  C(cl)  to  be  zero.  We  may  now  write  at 


x  =  b,  using  boundary  conditions  (b)  and  (c), 


3>(W  =  /?c<t.i  si  nViCr.b) «  ;B60  e 


-Yt 


(3.  15) 


^(b-o')-r  Ikb-t-o)  =  fl(a)  Vi  Cosh(Yib) 


-YL> 


(b+o)  •+"  ^  (o-to)  -  Y  ^ 


(3*  16) 


(3.  17) 


A  solution  for  C^(b)  will  yield  \(a_)  and  B(O-)  and  hence  a  solvation  for 
J(x.  G-)  and  a  formal  solution  for  (x.  z). 

..  Gc 

Best 


De*ine 
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D|_  =  ^.'(b+o)  — 


(3.  18  * 


.-1/2 


and  therefore  D  is  analytic  in  ‘7'>-Mm(k  ,  k  )  and  behaves  as  CJL  as 

*1*  w  4 

<X-*-oo  in  7“  >  -Min{k^,  k^).  Subtracting  equations  (3.  16)  and  (3.  17),  using 


the  boundary  condition  (d),  and  the  results  in  (3.  15)  gives 


0.  = 

where: 

L(c0  = _ S 1  n li  ( Y  jp) _ 

Y  Sin  In  (Y,  u)  YY|  C<?S  b(Yib)  (3.  20) 

The  function  L.'oL?  has  branch  points  at  k  and  -k,  the  branch  points  of 
Y  .  Choo  sing  the  branch  cuts  for  Y  as  was  done  in  section  2,  <2.  1 8 » ,  and 
shown  in  Fig.  2,  results  in  L(CC>  again  being  analytic  in  "k^^  '7"<C  k^.  This 
function  must  be  factored  into  a  product  L  L  .  The  function  L  is  analytic 

T  T 

in  7~> -k^  and  L  is  analytic  in  k^.  Again,  this  factorization  s  the  diffi¬ 

cult  step.  The  factorization  is  obtained  in  section  3  .  1.  3  by  a  limiting  proce¬ 
dure.  The  method  is  similar  to  the  one  used  in  section  2 ,  3  and  gives  the  fac¬ 
torization  in  a  form  convenient  for  numerical  work.  Multiplying  (3.  19)  by 
L4  (<3j  and  rearranging  yields 


ljb)  _  p1  i 
L«r!  'z  b  (a-ys,) 


D^.  Lw-)  -  E.W) 


lLk)-  +  Ew> 

L<d) 


n.  2ii 


where: 


Eio^-iT1  g  U-^L- 

bCoL-^i) 


(3.22) 


E(o-)  =  jJr^ 


b  (CL-yS,) 

E. .(aj  =-1Sn_i - [u®-)  -LCS,)] 

+  7*  bCt-^i  + 


(3.23) 


(3.24) 


Obviously  E  (a)  is  analytic  for  7"<  -k^  and  E  (d.)  is  analytic  for  7~> 

-Min(k  »  k  ).  Therefore  equation  (3.21)  holds  for  -Min(k^,  k  )  T «£. 

Min(k  ,  k  ).  The  left  side  of  (3.  21)  is  analytic  for  7“>  -Min(k  ,  k  )  and  the 
2  4  2  4 

right  side  is  analytic  for  /"<  Min(k^>  k^).  Therefore,  one  side  is  the  analy¬ 
tic  continuation  of  the  other  and  both  sides  of  (3.  21)  may  be  set  equal  to  J(CL), 
which  is  analytic  in  the  whole  QL -plane. 

-1/2  . 

In  section  3".  1.  3  we  will  find  that  L+  and  L  behave  as  d.  in 
'7">-k  and  '7"<C  k  »  respectively.  Therefore,  E  (  <3.)  behaves  as  CL  ^  for 

W  ta  “ 

k  and  E  (  cl  )  behaves  as  CL  for  7">  -Min(k  ,  k  ).  Using  these  re- 
suits  in  3.21  shows  that  J(<£)  behaves  as  CL  for  '7''  >•  -Min(k^,  k^)  and 
also  for  7~ <  Min(k  ,  k  ).  The  extended  form  of  Liouville's  theorem  proves 

A 

that  J(<X)  is  zero.  Hence,  we  obtain 


ACalJmUvlx)--/?  1  Ljfri SinUyjjO 

1  I  2  b  S,nVi  (Y,t) 


(3.25) 
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B  U;a)  =  B(cl)  e  =  -  t /?  z  >L--a)  --S 

*'  '*  b(OL-A.) 


Yb-Y* 


and  the  formal  solution  as 


&+Z  T  . 

a  — ^  CLZ. 


H,(x,a.)  e  da;  o<a<t, 


-  osttr 
co -fir 


cj)(x;2.)  = 


—00+27" 


with  -Min(k  »  k  X^T^Mintk  ,  k  )  in  (3.  27)  and  (3.  28). 

u  4  L* 


(3.  26) 


(3.27) 


(3.  28) 


3.  1.  2  Choice  of  a  Closed-Region  Structure 

The  chosen  C-R  structure  is  shown  in  Fig.  8.  The  formulation  of 
this  problem  is  identical  to  that  of  the  O-R  problem  in  section  3.  1.  1  except 
for  one  important  change.  The  boundary  condition  which  required  decaying 
waves  at  infinity  now  becomes 
f)  (£(a,  CL  )  =  0 

This  results  in  the  following  equation,  corresponding  to  (3.  14),  as  a  solution 
of  the  wave  equation  (3.  11). 


Jfr  =  g(a)  Sih k  [y (a-x^  +■  D(ac)  Cash  [y (a-x))  j  hi.  x  <  a 

Following  the  method  of  solution  as  outlined  in  section  3.  1.  1  gives 


(3.  29) 
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/?  e  K^Kfa)  S.'nU^i  *V 

'*  Mo.-/*,)  Sinh(r,b) 


(3.  30) 


D>  (*,a.)  ~  B (a) 5m (, \y (a-*)]  =  - ■/?  £  Kj^>)  K^l  Sihk[y 

'  %  t)  C<x-^,)  5i'rt^  Cyc) 


$(x.z)  = 


a+ZT 

r 


-fz^ 


—  ICLZ 

c^^a)  0  da: •  o<z c^t 


; 


c(>(x,z)- 


-eo+zT 

CO-t-iT 

n 


-]zrrr 


n 


-ZCLZ 


D,(xa l)  dd  , 


; 


-ic-t-zr 


(3.31) 


(3.  3Z) 


(3.  33) 


with  -Min(k  ,  k  )  zL/7'<  Min(k  ,  k  )  in  (3.  32)  and  (3.  33).  The  function  K  (a.) 

w  fr  L  t  T 

- 1/2 

is  analytic  for  7"  >  -Minfk^,  k^)  and  behaves  as  a.  with  cl-*»oo  for  '7"  > 

- 1  /2 

-Min(k  »  k  ).  K  (a.)  is  analytic  for  7"<.  Min(k  ,  k  )  and  behaves  as  a. 

M  T  “  U  1 

with  ai— ►CO  for  ^  <  Min(k  ,  k  ).  The  product  K  (  CL)  K  (<3-)  is  equal  to  K(&) 

u  *t  T  “ 

in  the  strip. 


K(cl)  = 


S'mUVjL)  (Vc^ 


Y  Smh(r,b)  Coj.Y(rc)  -t-Y,  Cask  (Y(b)«^UYc)  (3. 34) 


The  function  K(<£)  is  a  ratio  of  integral  functions.  It  is  actually  meromorphic 


with  poles  at 


= 


^TL 


(3.  35) 


with  1^  being  the  zeros  of  the  characteristic  equation  (transverse  ”  we  num¬ 


bers  of  the  waveguide)  for  the  inhomogeneously  filled  waveguide. 
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Co*Uc)  si*  Uti  +Cos(4D  S.J/c) 

A  A 

where: 

Note  that  the  zeros  approach  (nir/a)  as  1-^oc  . 

The  functions  K^(&-)  and  K  (d)  are  now  readily  obtained  by  using  the 
infinite  product  expansions  of  the  integral  functions. 


(3.36) 

(3.  37) 


f  = 


Sih  (kjk)  S\ hike) 


kS,'h(kdb)  Cos(kc)+  STih(kc) 

,  [h  /  _  ‘t  >  l/z 


n.  =  ( kL-  m 


Z\  Vi. 


«n  =  (k~- A 


Z  }  '/z. 
»  I 


(3.39) 

(3.40) 
(3.  40.  a) 
(3.  40.  b) 


(3.41) 


The 


X  (<£)  given  by  (3.41)  ensures  that  K^(OL  )  behaves  as 


-V2 

CL  as 


q_-»-  oo 


7~>  -Min(k^.  k^).  It  is  obtained  from  a  knowledge  of  the  asymptotic  form 


-'f- 


of  the  infin.*e  products  in  (3.  38)  by  means  of  equation  (2.  54.  a).  1  he  function 

K  (a)  is  equal  to  K^(-<J.).  This  function,  K(  a' .  and  its  factorization  will  yield 
the  factorization  of  L{<X),  equation  20),  via  limit  in  a  form  convenient  for 
numerical  calculations. 


3,  1,  3  Factorization  for  the  Open -Structure 

The  factorization  of  L)  QL ) »  equation  (3.  20),  will  be  done  by  a  limiting 
procedure  analogous  to  the  method  discussed  in  section  2„  3,  The  only  differ¬ 
ence  here  is  that  L(  CL  ’  is  more  difficult  ^compare  (3.  20  i  with  <?.  30  •  and  a 

closed  form  fur  *he  answer  is  not  obtainable.  However,  +he  form  o'  the  fac 


torization  obtained  readilv  .ends  itself  to  numerical  processing  and  hence  nu 
merical  resul's  for  +he  elec'romagnphc  Held  quantifies  of  inures*. 

The  relationship  between  *he  O-R  s’r-jc’ure  and  *he  C-R  s’ru<  tiro  .s 
seen  Ov  comparing  Fig.  '<  a;*1  8.  The  C  R  structure  .s  ob'ained  by  F.-**  ng  a, 
c:  — oo  while  maintaining  a  c  b.  Using  this  Iim.f  on  K(0L  ',  equa‘ion  <3.  34  • , 
for  any  a_  such  *haf  - k^<. k  ,  gives 


Li  tn  K(a)  =  Lihn 


^hl(Vib)  Sj.(Yc) 


a,c  — *•  <70 

a-c  =  b 


C“*’0°  )  Y’^rKv',  b)  Cos  h  (Yc  )  -f-  Yi  CoS  V  (Y,  t?)  ■S’lh  i\  LYc) 


y  j".  n  u  c  y;  i )  y,  Cc^  I,  Cr,  b ) 


=  L(a) 
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This  results  from  the  fact  that  for  any  CL  wi'hm  the  s’ rip  Y  has  a  posi’i’  e 


non  zero  real  part  Therefor*- 


(3.43) 


■>  t 


L-,w  =  i 

=~®  s,„Uyc)  (3-4i’ 

The  function  K(&)  can  be  expressed  in  a  convenient  (for  taking  limits) 
form  by  first  expressing  K^(C-)  and  K  (C£.)  in  a  form  like  ihe  one  used  in  sec¬ 
tion  2.  3.  In  particular,  we  have 

a.  b 

;2  r oo  ”  irrn 

■f0'7T(i  +  -f-)e  •  e 

Lr7=l  p*  J 

with  H  (CL)  defined  by  (2.  56),  (2.  57),  and  (2.  58).  Again  ^C(cl)  is  not  included 
ii  (3.44).  In  this  case  the  functions  F  (cj)  and  F  (u>)  are  the  charactt  ristic 

JL  L. 

equations  of  the  structure  shown  in  Fig.  8,  namely, 


1 

HCCd) 


(3.  44) 


with 


Fj(u)  =  Sin  (oic) 

f^(u)*4j£os('&Jc).5ih(aJ|k)  4-  OJfCosiu).  b)  rs(ooc-) 

’  .  2.  N,/2” 

CJ,  ”  (  ul  —  R.  > 


(3.45) 

(3.46) 

(3.  46.  a) 


The  contour  used  in  (2,  56)  must  enclose  the  proper  zeros  of  equa¬ 
tions  (3.45)  and  (3.  46),  that  is,  half  the  total  number  as  the  zeros  occur  in 
pairs  +  oj  .  Equation  (3.  46)  now  has  >he  possibility  of  zeros  that  give  rise  to 
surface  waves  and  again  zeros  whose  spacing  approaches  a  continuum  as  a,  c, 
approach  infinity.  This  is  clearly  demonstrated  by  considering  the  zeros  of 
(3.46)  when  the  loss  CTj  an3  ,  are  reduced  to  zero.  Under  this  condition 
there  are  two  possibilities,  real  roots  and  imaginary  roots.  The  imaginary 


roots  will  exist  if  is  greater  than  one,  which  is  the  case  being  considered 


here.  The  possibility  of  less  than  one  is  discussed  in  section  3.  3. 

The  real  zeros  of  equation  (3.46)  approach +mr /a  for  large,  real,  u 
and  also  the  spacing  between  the  zeros  becomes  infinite  jimally  small  as  a 
approaches  infinity.  For  imaginary  roots  letw  =  i  p  (p  a  real  number)  and 


equation  (3.  46)  becomes 

Fz(  if>)  =  Z  Sos  h  h  b)  -h  £oS  t)  S ih/i  ^ 

Op  V 

/C-l)^0  — j*  )  ^  *  Equation  (3.  47)  can  be  shown  only  to  have 

zeros  for  0  |  p  |<^  n/'/C  -  l'  k  .  The  number  of  zeros  is  discrete  and  there 

o 

may  actually  be  none.  How  many  real  zeros  equation  (3.47)  may  hive  is  de¬ 
termined  by  the  parameters  b,  ,  k^.  When  c  approaches  infinity  and  p  is 
positive  (3.47)  reduces  to 


Ay  n  ( ^ob)  *+*  'ja.  |  Coj 


(3.48) 


which  is  the  characteristic  equation  for  determining  the  surface  roots  of  a  di¬ 
electric  slab  backed  by  a  perfect  conductor  with  a  TE-polarization.  (See 
Collin  [i960],  p.  474).  Therefore,  the  imaginary  roots  of  (3.  46)  which  are  in 
the  upper  half  of  the  o-plane  go  into  the  surface  wave  roots  and  the  positive 
real  roots  go  into  the  continuous  eigenvalue  spectrum  as  a  and  c  approach  in¬ 
finity  for  the  open  structure  shown  in  Fig.  7.  h.’ctc  mat  equation  (3.  4b)  is 
slightly  different  from  equation  (3.  36).  However,  they  have  the  same  zeros 
because  1  =  0  or  1^  =  0  are  not  roots  of  (3.  3b).  A  zero  at  1  =  0  is  present 
only  when  the  parameters  b.  fc  ,  kQ  are  such  that  the  transition  point  of  a 
new  surface  wave  occurs.  We  will  pick  the  parameters  such  that  the  transi- 


tion  point  is  not  obtained.  When  the  loss  cf  and  &  are  reinsert'd  the  zer^s 
will  become  slightly  complex.  The  previous  real  roots  will  now  contain  a 
small  imaginary  component  and  the  imaginary  roots  will  contain  a  small  real 
component. 

The  contour  used  in  (3.44)  is  the  one  shown  in  Fig.  4  with  the  addition 

of  contours  to  pick  up  the  surface  wave  roots  in  the  upper  half-plane,  if  any. 

The  zeros  (ntr/a)  are  now  replaced  by  the  ~  ,  the  zeros  of  (3.46).  The  radius 

n 

r,  in  Fig.  4,  must  now  be  0  <1  r<ReL  ),  Im(-  )<  £  <£  k  .  where-  is  the 
°  o  o  2  o 

smallest  root  of  equation  (3.  46)  belonging  to  the  set  which  goes  into  a  contin¬ 
uous  spectrum  as  a  and  c  approach  infinity. 

That  K  ( Gt )  given  by  (3.44)  is  the  same  as  (3.  38)  can  be  verified  by 
using  the  calculus  of  residues.  The  only  singularities  within  the  contour  are 

l  • 

the  zeros  of  F^(cj)  and  F^(j).  F^  (u)  and  F,  (- )  do  not  have  any  poles  within 

X  • 

The  contour  integral  for  H  ( CL)  can  now  be  written  alcng  each  path 


and  gives  for  the  surface  waves 


60 


with  H  (d  )  given  by  (2.  o5).  The  function  K  (Ct)  -  K.  CL).  Hence, 


1 
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r  _  cl 

l/j. 
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03  .  “in7r 

TT(i+-f-)e 
_"=i  pn  J 
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_T»-I  1k+pr  J 

hV) 
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•  T  • 

'0 


"co 

TTO-f- 
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L-1  w^] 

j 

.  e 


Taking  the  limit,  of  (3.  51',  and  using  (3.  42),  (3.  39)  yields 


L«  hr\ 

a.jC-*-  cd 

=  b 


KCol) 


(3.  51; 


P  ab  i 

ikl 

Sin(kib)  * 

oo  ^  t  h  TT 

ttci+4-  e 

n*l  rn 

■i  1  —  - 

• 

TT(i-fn)e^ 

Lr  =  i  Pn  J 

-£k  -Sth  +■  Co^(  kdt>) 


h  •  |  l/k+f3* 


M 

TT(i- 

fl  =  ‘V 


H(cO  h-HC-cL) 


<  3.  52  i 


where  p  is  now  a  surface  wave  root  (positive  roo1'  of  (3,48)  .  M  +he  num- 


n 


ber  of  rootsj 


HCgl)  =  Li*!  HCO 


tx,c  -+ 

9C-  C  c  b 


(3.  53) 


with  H  (CL)  given  by  (2.  6C). 


Define  (  £  >  0< 
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h-  -f.  Gr(f  =  Lim  311  LJ. 

ZT  *  o^jC+cs  Z7r  i 

o.-C  *  b 


-L-  L;m  JVil^Ll) 


&H«) 

Ft(r-u) 


r 

j 


r 

j  s;„(w,l>)  +  ub 

jj*  Cos(w,b)+e5h(^t|-.  i 

CosCcJ.b) 

j  zri  J 

__W  Sin  Cw,b')  -  i  U,CoS  CtJ,b)  j 

! 

i 

and 


(3.  54) 


i—i  _ _ 

cl, C  — *“  cC  --7T  £ 

-  Car  b 


r 


1|r)(^.g  T  Cs^0'^  "  2^in(^k)|+  Cos(y.b') 


-^2  rio5in(tJ(b) -{"  t  OJ,  CoS  (w,  b)^] 


(3.  55) 


so 


w-^+U 


<x> 


CJ-l€ 


<x>+?€ 


cc  r 

H(cO  =  -^  4?(a)w)cJrj  +1  fca/^G^u;)  dw  —  f  (<Z}lo)Gz(u)  doj 

Wso  COxC-id  .\  ~  >%£. 


Cd  xO+Z€ 


(3.  56) 


with  Imu  <  6<  kv  In  this  form  the  terms  due  to  G.  and  G?  can  be  consider- 
ed  the  perturbation,  due  to  K  varying  from  1.  The  equations  (3.  52)  and 
{  '  ”6)  give  L(OL)  in  a  suitable  form  for  obtaining  L^(OL)  and  L  (OL). 

The  choice  lor  L^(CL)  may  be  made  as 
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'/2 

_£k| 

Sin(KiQ 

9 

roo 

TT(i+|-)e 

M  =  l  pr  j 

"*tis  SirtCfej  W (l^b) 

(3.  57} 


L_(a)  -  L  +(-d)  (3.58) 

The  function  L  ( cX )  is  analytic  for  'T'  >  -Min(k_f  kj.  This  results  from  the 
+  L  4 

fact  that  £  ,  eouation  (2.40),  has  a  positive  imaginary  part  great  er  than  k  , 
n  4 

and  therefore  the  infinite  product  is  analytic  for  7~  >  -k^.  The  function  H(CL  ), 
equation  (3.  56),  is  analytic  for  as  f(  (t ,  w)  and  G^(gj)  and  G^(w)  satisfy 

the  conditions  of  theorem  A,  page  11,  Noble  [1958]  and  restated  in  the  Appen¬ 
dix. 


The  solution  of  the  Wiener -Hcpf  type  equation,  (3.  21),  requires  that 

L  ( <X  )  and  hence  L  (rs.)  have  algebraic  behavior  for  large  CL  .  Therefore, 

-  )C  ( C2 ) 

the  term  e  must  be  included  in  (s.  57)  to  ensure  this  behavior.  The 

asymptotic  form  of  (3.  57)  will  dictate  the  choice  of  X  (CC). 

The  infinite  product  behaves  as 


(3.  59) 
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for  "7“>-Min(k2»  k  ).  The  finite  product  behaves  as 


M  ,  cL 

TT(l  +  jfr=r'e 

»=<  flf+p; 


M  OL 

g^ww 


Also 


<30 


~7T(\  ^  ^  ~  2  —  2  ^  In  /^S 


W=o 


which  is  obtained  from  the  work  in  section  2.  3. 


This  leaves  the  term 


(3.  60) 


(3.  bl) 


CD  -t-£ 

r 


<30  f2€ 

-fc^jCj)  G{  (u)  <Juj  -  J  -pCcLjCj)  GJd)  dot 

(i)zo+i6  (3.62) 

Consider  it  in  two  parts  by  referring  to  the  definition  of  f(a,  u)  given  by  (2.  57). 
One  term  is 


r 


ff>  -  26: 


<35+ 


Ha)  -  -  aJ  ( kz-  w4) ' ^  Gj(^)  J  cj  -f  ( fe*-  Gj/w)  Jw 


<x>=o  —26 


Wro+l6 


(3.  63) 


v/hich  behaves  as  ((X)  .  (constant).  These  integrals  converge  because  Gj(u>) 
and  G^co)  go  to  zero  for  large  u>  on  the  contours  indicated  in  (3.  63).  The  re¬ 
maining  terms  of  (3.  62)  are 


03-2.6 
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lid) 


CD  +  tS 
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*  ('  ~  ] lh  ('  +  ffep)  ^  (3‘ 641 


w=0— 2>^ 


,0i)ao  +2  6 
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Changing  the  variable  in  the  second  integral  of  (3.  64)  to  -u  and  using  the  def¬ 
initions  of  G^(gj)  and  G^w)  given  by  equations  (3.  54)  and  (3.  55),  respectively, 
gives  for  I ^  ( dt  ) 


(3.65) 


The  path  of  integration  may  be  closed  in  the  lower  half  plane  since  G^(w)  goes 
to  zero  on  the  infinite  semicircle.  The  integrand  of  (3.  65),  however,  may 
have  poles  at  co  =  ip  where  p  is  the  surface  wave  pole  at  equation  (3.48),  and 
also  the  branch  of  nT k  +  go  is  in  the  lower  half  plane.  Therefore, 


The  branch  line  integral  of  (3.  66)  is  equal  to 


G'j(u>')  dw 


(3.  66) 


-ico 


Cofistiht 


-Ip 

as  oo  .  Therefore 

X.Col) 


h  (  <XM) 


(3.67) 


(3.  68) 


for  7~  >  -Min(k  »  k  ),  Combining  these  results  gives 
c$  *± 
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IM)  ~  <* 2  rr  +  lh(kt)l+  &  +  Tea.)  +£  ,  I1—  ti 

r  |  77-  L  VzTT/J  Z  nsi^+Ptt 


(3.69) 


Hence,  the  choice  of  X  (CL)  as 


*  w  =  ExP  [p  +  ImS)] +  ^  +  IC<L)  +  al 


(3.  70) 


will  result  in  L+(C7.  )  remaining  analytic  for  '7"  >  k^)  but  now  behav- 

-1/2 

ing  algebraically,  namely,  OC  ,  as  (L-+-00  in  this  region.  The  function 
1j  (CL)  is  given  by  L+(-a)  and  hence  L  (<2)  behaves  as  CL  as  CL-*-5’'*1  for 
7“<  Min(k7»  k  ).  The  product  L  (d  )  L  (<£.)  is  still  L(a. )  in  the  strip  because 
of  the  fact  tnat  X  (~  OL)  -  ~X  (3.)  (refer  to  the  definition  of  X(QL)  given  by 
equation  (3.  70)  ). 

Reducing  losses,  O',  and  £  ,  to  zero  once  L  (cL)  is  obtained,  yields 

X  Cm  i 

a  simplification  in  the  equations.  Setting  <3*.  and  to  zero  results  in  k  =  k  , 

k  =  n/^"1  k  .  Also  in  equations  (3.56)  and  (3.63).  Now  G  (u>)  =  -G  (w) 

do  2  1 

with  w  real  (refer  to  equations  (3.  54)  and  (3.  55)  defining  G.(w)  and  G  (co),  re- 

X  Cm 

spectivcly)  and  the  loss  reduced  to  zero.  Therefore,  we  obtain  for  (3.  56)  and 


(3.  63),  under  the  condition  that  the  loss  is  zero, 


HCd)  *r  J  -f  -h  CricJ)j  d< 


(3.71) 


—  Vz 

XU)  -  -  (k*-wz)  C*  (<*>)  dw 

c 


(3.  72) 
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where: 


G(<d)  =  Z  Re  &Aui)  =  k. 

•  -rr 
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f  b  —CoS  (<*•'; t>)  S'i n  (w,b) 


OJ2-  ^ 
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LW,  I 
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'.r  [co2-  Sin  (u,  b)  +  6dj*  Co/VoJ,b)  J 

/  i  2-  \ 

(  —  6o^” -f—  Cia-\)  k0  / 


(3.  73) 


(3.  73.  a) 


and  5  .  is  the  contour  shown  in  Fig.  6.  Convergence  of  the  integral  in  the  de- 
finition  of  I(  OL ) »  (3.  72),  is  ensured  as  G(co),  (3.  73),  goes  to  zero  as  cj->oo  . 

A  form  of  L^(cl)  that  is  extremely  convenient  for  numerical  work 


when  the  loss  is  reduced  to  zero  is  given  by 
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l_(al 
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TT(^+a)e 
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Vz 


Hte)-X(a) 
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(3.  74) 


H((aO=l  |  ^ 

**  A  * 


(3.  75) 


with  X  (<X)  defined  by  equation  (3.  70)  with  k  =  k  ,  I(CL)  by  equation  (3.  72), 

G(oj)  by  equation  (3.  73),  p  by  equation  (3.40),  V  by  equation  (3.  12),  p  the 

n  i  n 

real  positive  zeros  of  (3.48),  and  by  setting  k  k  ,  k  =  k  wherever 

d  o  o 

they  occur. 
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3.  1.  4  Solution  of  the  Problem: 

The  evaluation  of  the  field  quantities  of  interest  is  obtained  by  the 
Fourier  inversion  integrals  given  in  (3.  27)  and  (3.  28).  The  function  L+(ct) 
is  given  by  (3.  74)  and  L  ( CL )  =  L+(-cL).  Within  the  waveguide  section  of  the 
structure,  z  >0  and  0  <  x  <  b,  the  modes  are  obtained  from  equation  (3.  27). 
The  contour  may  be  closed  in  the  lower  half  plane  enclosing  only  pole-type 
singularities.  The  calculus  of  residues  gives  the  modes  directly,  in  particu¬ 
lar,  the  lowest  order  reflected  mode  (only  reflected  mode  carrying  average 
real  power)  is 


nr 


/» 


Hence,  the  reflection  coefficient  is 


(3.  76) 


_ - 3f—  L  W 

Zb  ft  (3.77) 

The  far  field  is  again  obtained  bv  the  use  of  the  Huygen  equivalent  source  in 
the  aperture  x  =  b  +  0  in  a  manner  analogous  to  that  described  in  section  ?.  4. 
The  only  difference  is  now  we  have  an  equivalent  magnetic  current  in  the  aper¬ 
ture.  The  results  are 


s,w«). 

k,  C os  6  -h  8\ 


(3.  78) 
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I  2  2 

with  p  =  -|x  +  z  ,  0  «£  ir,  and  kQ  =  uW  jj  q  6 

In  the  results  given  in  (3.  77)  and  (3.  78)  the  loss  has  been  made  negli¬ 
gible!  in  fact*  zero.  The  loss  set  at  zero  results  in  k  =  k^  and  k^  =  \ f/0 k^  a  id 

hence  equation  (3.  78).  The  loss  will  also  be  set  at  zero  in  the  results  that 

follow. 

The  average  real  power  reflected  in  the  waveguide  and  the  average 


real  power  radiated  in  the  space  wave  are 

Normalized  Reflected  Power  =  W  ,  -  IrI 

ref 

Normalized  Radiated  Power  =  W  ,  = 


(3.  79) 


z  nr 


4* 

k0Cos  p 


(3.80) 


The  normalization  consists  in  the  incident  power  being  set  to  one.  Recall  that 
the  parameters  are  chosen  so  that  only  the  lowest  order  mode  propagates  in 
the  waveguide.  Therefore,  the  higher  order  reflected  modes  carry  zero  aver¬ 
age  real  power. 

The  structure  of  Fig.  7  has  the  possibility  of  the  existence  of  surface 
waves.  The  surface  wave  modes  for  b  ^  x  are  obtained  from 


,  f  rW  -  Vx-?aZ 

[u,z)=.  jjygj  s e _ jg--  z<o 

i  }  2-t  J  (a.-*,)  >  .  .  .nU(v'.k)  +  y,  Gwfcfak)]  Ua:) 

-a-rtr  r  r  (3.81) 

Use  of  L.  (  CL)  -  h(  a.)/L+(  ql)  is  made  in  (3.  8  i ).  The  surface  waves  are  given 
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by  a  j  the  positive  zeros  of 
n 

Y  Sir'll  (y.b)  •fV.  (2osV.  (Yt  L>)  (3. 81. a) 

with  c X  real  and  k  <  &.  k  .  The  possibility  of  surface  waves  exists 

n  o  n  o 

because  y£  >  1  is  being  considered  here.  The  residues  of  the  integrand  of 
(3.  81)  at  the  surface  wave  poles  give  the  surface  wave  modes  as 
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" <S,Z 


2.  Resn  *  e  •  b<*  Z<o 

y  b  r  =  i  >  J 


(3.82) 


where: 
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esh  — 


^  pn 


0,nb) -  a^CisC^t)  Sih(^y  -V-  Wah  p„  L 


(3.83) 


p  (0  <  p  <  \I/C  -  l'k  )  is  a  positive  root  of  (3.48)  with  the  loss  set  to  zero, 
n  n  o 

M  =  number  o.  surface  waves  (number  of  positive  real  zeros  of  (3.48)  ) 


j  z  n 

CL  =  Ik  +  p 
n  f  o  n 

=  V feo  P n  > 


(3.  84) 


(3.  85) 


The  surface  wave  modes  for  0  <  x  <  b  are  obtained  from 


e  -iM'’ 

j  *t» 
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Sink  (Y,x) 


-2U.Z 


UaXa-^)  [Y  5;„>,[Y,b)  -t-  Y,<T:,k  «t)] 


b  z< 


_«  +  1 7"  '  “*  (3.  86) 

The  residue  of  the  integrand  of  (3.  86)  at  the  surface  wave  poles  again  gives 


for  the  surface  wave  modes. 
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‘'04SC1b  Z<°  (3.87) 


Sf  h  (V)mL)  ^ 

The  surface  wave  modes  are  orthogonal  in  the  following  sense: 
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(3.  89) 


Hence,  the  total  real  average  power  carried  by  the  surface  wave  modes  is  the 
sum  of  the  power  in  each  mode.  This  gives  the  normalized  real  average  pow 
carried  in  the  surface  waves  as 


er 
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An  inspection  of  the  results  given  in  this  section  shows  that  once 
L+(tf)  is  known  (5  real),  the  radiation  pattern  of  the  space  wave,  the  power 
reflected  in  the  waveguide,  the  power  carried  in  tne  surface  waves,  and  the 
power  carried  in  the  space  wave  can  be  determined.  These  numerical  calcu- 
-ations  are  discussed  and  the  results  presented  in  section  4. 

The  radiation  of  a  truncated  paraJlel  plate  waveguide  in  free  space. 
Fig.  1,  with  an  E  in  ident  may  be  obtained  i  om  the  results  of  section  3.  1  bv 
setting  the  relative  lelectric  constant  to  ore.  ThiS  value  of  .he  dielectric 


constant  removes  any  surface  wave  phenomenon  and  hence  any  equations  foi 
surface  wave  quantities  are  set  to  zero. 


3.  2  TEM  Excitation  of  a  Surface  Wave  Structure 


In  this  section  we  will  obtain  the  excitation  of  the  surface  wave  struc¬ 
ture,  shown  in  Fig.  7,  for  a  TEM  mode  incident  in  the  waveguide.  The  inci¬ 


dent  field  is 


H  =e  d  Z;  0  <  x  ^  b,  V  Z 


with  k  given  by  equation  (3.4).  The  formulation  of  the  problem  follows  ex- 

d 

actly  the  one  given  in  section  2  except  that  now  tha  presence  of  the  dielectric 


must  be  taken  into  account.  The  results  of  the  analysis  are 
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where. 
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(3.  97) 


B(.,o=_i!=ttWLu) « 
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'fwP  Y  (a  -  kj) 


(-5.  9- ) 


c/ith  /Z  }  s  Y^/Y^  and  V1  defined  by  equation  (3.  12)  “_nd  Y  by  (2.  lo). 

The  function  L(  CL )  for  this  polarizatior  of  the  incident  field  is 


.(cl)  = 


vy,  ^Uy;k) 


y.  Sihk(Y,b)-bY/C,  CosUyb)  (3.99) 

This  function  is  factored  by  th^  rr ethod  of  section  3.  1.  3  by  using  the  related 
closed- region  function 


with  the  result 


_ YY,  Si  h  ii  (  Yc  )  sSi'h  b  (fib)  _ _ 

f,  •S  *  h  Cy,  b)  Cos  )n  (Yc)  +-KY  Cos  Ii  Yb)  -Sin  h  CYc) 


'/r* 

L^)  =  (cL+fe)"  (.CL+ki)  MjflC) 

-f-  1  T” 


(3.  100) 


(3.  101) 
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i*JJL(,+;s^e  J*e 


H(tt')  —  X(rt> 


[  fej  S;n(kjb)  -zK.k  C«‘^k)].lj*  (l+ jj^—j)  e' 


(3.  102) 


L(<0  =  ( <x-k)Ccx— k.)  M(-a) 


(3.  ?03.  a) 


and  k  is  defii.-d  by  equation  (2.  4)  and  p  by  (3.  40).  The  numbers  p  are  the 

n 

surface  wave  zeros  given  by  the  zeros  of  the  characteristic  equation. 

/£.  o  Cas(kf>ij  =  p(  (3.103) 

R  =  -kZ~f>Z)'/*'  (3.104) 

The  function  Hi  a)  is  still  given  by  (3.  56),  X  (  CL)  by  (3.  70),  and  1(0.)  by 

(3.  b3),  but  now  the  functions  G  (w)  and  G7(co)  become  for  this  problem 

r  .  x  ,  -2/f ,  Cts(w,t>)  "t-  73”  +(o\?  \c*s(v,b)  d-i/k(Tj-  3ih(^'b)J 

Lriti))  =  £_  — - - -  -i - 1 - (3#  iQ 

Z7T  2,7T:  [-la,#' +  0 j.^inCm,)"] 

£t>,  "  (  W  H-  fe  /  (3#  lo6 

-  —  6^(-or)  (3.107 

This  results  in  b^(&,  being  analytic  and  benaving  as  as  Q  ,-lqq  for 

^  >  -Mintk^,  k^,).  Likewise  L  (cl)  is  analytic  and  behaves  as  as 

Or*<»for  r"<  Min(k_,  k  ), 

2  H 

The  form  of  1^  a. )  that  is  convenient  for  numerical  calculations  when 
the  loss  o'  ^  and  are  reduced  to  zero  is  given  by 


(3.  105) 


(3.  106) 
(3.  107) 
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Y,  Slh  Vi  (Y,k)  ffcY  L?o$ll  (y,  t)j 
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(3.  108) 


with  ^  the  contour  shown  in  Fig.  6  and  G(co)  given  by 

/C  b)  CosCtJ.  b)  jjoJ7  ""  ^j]  £  w1  Is 

'TT  '\K.7‘  q7-  Cos  (CJ,  b)  +  Olf 


(3.  109) 


GYco)  — ~  — 


,  .  2.  \ 

61,  -  (  a/>  c/C- 1)  kc  ) 


(3.  110) 


(3.  ’ll) 


The  )(  (  OL)  is  given  by  (3.  70)  with  k^  -  YJ*?  k^,  k  =  k^  replaced  throughout. 
The  I(  OL)  used  in  )C  ( <*•)  now  becomes 


XCa)  =  —  ^  Ck^-  * J1")  Grt<jj)  doo 


(3.  112) 


with  G(<o)  given  in  (3.  110).  This  integral  converges  as  G(u)  behaves  as  Sin(w) 

for  large  co.  The  p  become  the  positive  real  zeros  of  (3.  103)  with  k  =  M00  k  , 

n  do 

k  =  k  ,  fc.-K- 

o  1 

The  results  of  interest  for  this  problem  are  obtained  as  was  done  in 
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section  3.  1.4  and  are  merely  stated  here.  The  normalized  (incident  energy 
set  at  one)  reflected  power 


(5.  113) 


Th;r  follows  since  the  parameters  are  restricted  to  values  such  that  only  the 

1/2  . 

lowest  order  mode  can  propagate  in  the  waveguide.  That  is,  /C  k  ^  ir/b. 
The  normalized  power  radiated  in  the  space  wave  is 


w  ,  =  £ 

^  £5rEIT 


L/\>| 

ii 

ii 

e-o  I 

l_CK,Cose) 


with  the  radiation  pattern  given  by  (0  ^^tt) 


1  Hj  * 

LCfcs  Cos  ©) 

* p  y 1 

k0  Cos  e  +  K k 

The  normalized  energy  in  the  surface  waves  becomes 


(3.  114) 


(3.  115) 


where  M  is  the  number  of  surface  waves  and 


Dresn 


.  p> _ 

(S-K>k„)-  U>)  •  D 


(3.  117) 
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D-a„pn  Si„(^)+kh[Jhkc«^W+a»klc^Al=)+l>ahp^/c^''^»l»)  (3  118) 

a„=(k.+F>J  <3-119) 

k-  (ot-i) feT  —  p*)'/*'  a<  <3-i2°) 

These  numerical  calculations  are  discussed  and  the  results  presented  in  sec¬ 
tion  4  along  with  those  for  the  TE-polarization. 

3.  3  Excitation  of  an  Incompressible,  Isotropic,  Plasma  Slab 

The  structure  in  question  is  still  given  by  Fig.  7  but  with  one  change. 
The  dielectric  between  0  <  x  <  b  is  replaced  by  an  incompressible,  isotropic, 
plcisma  medium.  This  plasma  medium  behaves  as  a  dielectric  with  a  relative 
dielectric  constant  less  than  1.  This  results  from  the  fact  that  the  relative 
electric  constant  for  this  medium  is 

£  “  vl  -X)  (3.  121) 

X  (3-  122> 

with  the  plasma  frequency  and  u  the  wave  frequency,  for  example  Budden 
[1964]. 

Values  of  the  relative  dielectric  constant  less  than  zero  are  not  of  any 
interest  as  it  is  impossible  to  have  propagating  waves  in  the  waveguide.  How¬ 
ever,  there  are  propagating  modes  for  the  relative  dielectric  constant  between 
zero  and  one.  It  can  be  shown,  for  this  range  of  fc  ,  that  the  strurture  will 
not  support  surface  waves. 
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The  analysis  of  this  problem  may  be  extracted  from  sections  3.  1  and 
3.  2  by  considering  the  effect  of  0  <  /C  1  on  the  analysis.  It  is  found  that  the 
only  change  occurs  in  G(u),  given  by  (3.  73)  for  the  TE  case  and  by  (3.  110)  for 


the  TEM  case.  Recall  thatco^  in  these  equations  is 


“i=  (k-'i^  +  u2)172 


(3.  123) 


with  u  varying  from  zero  to  infinity.  The  quantity  (/C-l)  k  is  now  negative. 

o 

Hence,  for  values  of  w  <  (1  -  /clkQ>  becomes  imaginary.  The  choice  oi  the 

sign  of  the  imaginary  number  is  immaterial  as  G(u)  is  not  a  multivalued  func¬ 
tion  of  w.  Therefore,  the  results  of  sections  3.  1  and  3.  2  stand  unaltered  for 


the  case  of  0</C<l.  However,  all  surface  wave  phenomena  is  non-existent  for 
this  range  of  K  and  the  equations  in  those  sections  must  be  interpreted  accord 
ingly. 


3.  4  Discussion  of  the  Method  for  the  D  electric  Slab  Structure 

Basically  the  related  C-R  structure  was  used  only  to  obtain  the  O-R 
factorization.  However,  the  O-R  and  C-R  solutions  may  be  related  as  was 
done  in  section  2.  5  in  the  following  way.  To  be  specific,  we  will  discuss  the 
TE  case. 

The  function  L+(  CL)  given  by  (3.  57)  is  the  limit,  a,  c-**0«  while  a-c  =  b, 
of  K+(  CL)  given  by  (3.  38)  if  we  include  the  )C  (<L),  (3.  41),  in  K^(  cL)  (recaJI  *nat 
makes  K+(&)  behave  algebraically).  This  is  true  if  the  limit  of  (3.41)  is  (3.  70). 


W e  may  write  (3.41)  as 
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^a)=^Co)+a£[(ki-C)/l  - 


(3.  124) 


Convergence  of  the  series  is  ensured  because  1  =  mr/a  +  0^1 /n)  for  large  n. 

The  function  Q((L)  is  equal  to  ( GL )  in.  section  2.  5  and  is  given  by  (2.  104). 
Recall  that  the  1  are  the  roots  ol  equation  (3.  36).  Taking  the  limit  of  X  (  ^) 


given  by  (3.  124)  yields 


Li  m  4-  2  ^  Wisk.)  4-  Sk  +3L(<l)  *H  Ct  2.  p2f  (3.125) 

ol.c-^oo7  77-  7T  \%nr!  2,  rh  v 

a-csb 

with  I(  OL)  given  by  (3.  63),  The  limit  of  Q((£)  had  already  been  obtained  in  sec¬ 
tion  2.  5;  refer  to  equation  (2.  109).  Now  equation  (3.  125)  Is  JC  (®r-)»  (3.  70)  for 
the  O-K  problem.  Therefore,  the  limit  of  K  (&)  is  L+(<L).  Hence,  one  may 
now  show  th*t  cfte  limit  of  C^x,  CL  )  and  Dj(x,  OL),  given  by  equations  (3.  30) 
and  (3.31),  respectively,  for  the  C-R  structure,  become  A^(x,  GL  )  and  B ^ (x,  GL), 
given  by  equations  (3.  25)  and  (3.  26),  respectively,  for  the  O-R  structure. 
Therefore,  the  complete  C-R  solution  becomes  in  the  limit,  the  O-R  solution. 
This  means  that  instead  of  formulating  the  O-R  problem  one  could  obtain  the 
solution  by  formulating  a  related  C-R  problem  and  taking  the  limit  ol  its  solu¬ 
tion  as  dictated  by  the  physics.  This  phenomena  is  what  Talanov  [  1 9 5 9 v]  sug¬ 
gested  would  happen.  In  fact,  he  solved  the  C-R  problem  for  a  TEM  exci¬ 
tation  and  calculated  the  energy  in  the  slow  waves  of  the  inhomogeneously 
filled  guide  and  claimed  that  their  value  as  a  — ►  <X>  is  the  energy  in  the  sur- 
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face  waves. 

In  section  2  we  obtained  the  factoriz  tion  for  the  C-R  problem  in  a 
particular  form.  The  advantage  of  doing  this  should  now  be  obvious.  As  a 
and  c  become  larr  e,  the  series  which  icr  ir red  could  be  recast  as  integrals 
providing  the  spacing  between  adjacent  zeros  of  the  characteristic  equations, 

F  (oj )  and  F  (u),  was  known.  This  representation  yielded  the  information  as 

1  L 

the  function  obtained,  -b  +  G(gj),  by  taking  the  limit  of  the  integral  representa- 

ir 

tion,  is  the  desired  information.  In  particular,  -b  +  G(o),  should  be  the  differ- 

ir 

ence  of  the  inverse  of  the  spacing  of  the  zeros  of  the  characteristic  equations. 
That  this  is  so  is  verified  in  section  4  witn  numerical  calculations. 
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4.  NUMERICAL  RESULTS 


The  usual  field  quantities  of  interest  ir  problems  of  this  type  are  ob¬ 
tainable  once  L+(d)  is  known  (drcal).  That  is,  the  average  power  radiated 
in  the  space  wave,  reflected  in  the  waveguide,  and  trapped  in  the  surface 
waves  (if  any)  and  the  radiation  pattern  of  the  space  wave  are  given  as  func  • 
tions  of  jL^(tf)j  .  This  can  be  seen  by  referring  to  equations  (3.  79),  (3,  30), 
(3,89),  and  (3.  78)  for  the  TE  polarization  and  (3.  1 13),  (3.114),  (3.116),  and 
(3.  115)  for  the  TEM  excitation. 

The  function  |  L+(tf)|  is  given  in  convenient  form  for  numerical  pro¬ 
cessing  by  (3.  74)  and  (3.  108)  for  the  TE  and  TEM  polarizations  respectively. 
To  indicate  what  is  involved  in  the  evaluation  of  j  L^td)]  ,  a  close  IcoK  at  the 
TE  cise  is  made.  The  evaluation  of  |L+(<f  )|  for  the  TEM  polarization  will  be 
similar.  Taking  the  absolute  value  of  (3.  74)  for  real  arguments  gives 


1  LW)| « 

Vk 
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3,(^-s) 
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The  real  pa»-t  of  is  given  by 
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Pell  lijgZ+g ) 


i  4-  6*(w)  ]  d  a) 

V  7T 


Ih  k-^l£i  >  •(-);_ -t6r(ui)i  du)  +&h  +l?el^) 


r*"“ -ffU 

This  result  s  obtained  sinre  the  square  root  term  is  purely  imaginary 


w  >  k0‘  The  definition  of  0  ,  (3.  40),  and  the  fact  that 


we  are  considering  the 


case 


where  only  one  mode  propagates  in  the  waveguide  ^  } 


gives 


It,  (/*" ) 

ft +6 

r-*' 

Therefore,  equation  (4.  1)  becomes 


(4.3) 


(4.  4) 


with  given  bv 
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HC«)  =  l 


(m-'O 


°'u' 


(4.5) 


and  G(cj)  <5iven  by  (3.  73).  The  convenience  of  L+(Cf4  as  given  in  (4.4)  is  low 
obvious.  We  are  left  with  only  a  finite  integral  to  evaluate. 

./hen  |^j  is  less  than  k  ,  the  integrand  of  equation  (4.  5)  becomes  sin¬ 
gular  at  <j  =  ^ This  results  from  the  fact  that  we  have  reduced 

the  loss  to  zero.  However,  this  is  an  integrable  singularity,  as  one  expects, 
and  may  be  handled  in  the  following  way.  Rewriting  equation  (4.  5)  for  0<(££  k 
gives 
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(4.6) 


The  first  integral  m  equation  (4,  6)  is  no  longer  singular  at  so  and  is  conven¬ 
iently  handled  by  a  digital  computer.  The  second  integral  in  {4.  6)  can  be  ob¬ 
tained  in  closed  form.  For  negative  values  of  d  we  have,  H  (6  )  -  -H  (Id*  ). 

z  z 

Therefore,  H^(d)  is  obtained  from  (4.  6>  for  all  d  and  hence  [  L+(d  )|  is  conven¬ 
iently  calculated.  The  characteristic  response  of  the  structures  for  both  TE 
and  TEM  excitations  will  be  giwn  in  graphical  form  in  what  follows. 

In  going  from  a  series  to  an  integral,  which  occurs  when  obtaining 
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L^(q.)  and  L  (OL)  by  letting  a  and  c  approach  infinity  in  the  representation  of 
K(QL),  we  needed  knowledge  of  the  spacing  of  the  zeros  of  the  characteristic  e- 
quatic.is.  We  concluded  that  the  function  b/ir  -  G(~)  is  the  difference  of  the  in¬ 
verse  of  the  spacing  of  the  zeros  of  the  characteristic  equations  F^(cj)  and 
■'7(cj)  as  a,  c  — +oo  .  Tor  example,  the  characteristic  equations  are  given  by 
(3.45)  and  (3.46)  and  G(co)  by  (3.  7 3)  for  the  TE  polarization.  This  conclusion 
r  tv  te  verified  by  actually  finding  the  zeros  of  F  (cj)  and  F  (t~)  and  seeing  if 

X  Cm 

indeed  the  zero*  behave,  as  a  and  c  approach  infinity,  in  a  manner  given  by  the 

function  b  -  G(uj).  The  result  oi  this  is  shown  in  Fig.  9  for  the  TE  excitation. 

■sr 

One  can  see  that,  for  e/b  equal  to  80,  the  curve  obtained  from  a  knowledge  of 
the  zeros  is  identical  to  the  theoretical  curve  given  by  b/V-G(u,).  For  c/b  -  8 
there  is  some  difference,  as  expected.  That  is,  only  in  the  limit  as  a—1 ►oo 
does  G(u)  hold. 

The  response  of  the  structure  to  a  TE  polarized  source  is  shown  in 
Figs.  10  to  15.  The  relative  dielectric  constant  ( /C )  has  three  district  ranges: 
a)  the  plasma  phenomena  with  0  <  £  <  1.  0;  b)  free  space  radiation  with 
/C  —  1.0;  c)  the  surface  wave  phenomena  with  K  >  1.0.  The  response  of 
the  structure  fer  "alues  of  /C  in  each  of  these  ranges  is  given  in  Figs.  lO,  ii, 
12,  and  13.  Recall  that  these  are  normalized  values  with  the  incident  energy 
set  at  1  or  100  per  cent  and  the  maximum  far  field  value  set  at  one. 

An  inspection  of  Figs.  10  and  11  shows  that  it  is  possible  to  have  all 
the  incident  power  radiated  in  the  space  wave.  This  ’s  also  true  for  the  sur  • 
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Fig.  9  Comparison  of  b/tr -G(i~)  given  by  (3.  73)  with  the  results  obtained  from  a  knowledge  of 
the  zeros  of  (3.  4  5)  and  (3.  46). 
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Power  distribution  and  far  field  patterns  for  an  isotropic,  incom¬ 
pressible,  TE  excited,  plasma  slab. 
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Fig.  13  Power  distribution  for  a  TE  excited 
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face  wave  structure,  shown  in  Fig,  1?,  until  the  structure  starts  to  suppor- 
surface  waves.  Once  surface  waves  can  be  supported,  a  transfer  of  power 
from  the  space  wave  tc  the  surface  wave  occurs.  The  structure  will  suppor* 

surface  waves  only  when  \"/C - 1'  k  b>tr/;>.  The  incident  power  may  all  be 

o 

transferred  to  the  surface  waves  by  a  suitable  choice  of  parameters  as  seen 
in  Fig.  13.  There  are  no  far  field  patterns  plotted  in  Ftg.  13  as  the  power  ra¬ 
diated  in  this  case  is  negligible  compared  to  the  power  tamed  in  the  surface 
waves.  Fig.  !3  also  shows  the  effect  of  two  surface  waves  existing  on  the 
structure.  Thi"  is  the  maximum  number  that  the  structure  can  support  since 
we  assumed  that  the  parameters  are  such  that  only  the  lowest,  order  mode  can 
propagate  in  the  waveguide.  The  effect  of  holding  k^b  constant  and  varying  1C 
is  shown  in  Figs.  14  and  15. 

The  r  aponse  of  the  structure  to  a  TEM  polarized  source  is  grvtn  in 
Figs.  16  to  19.  The  trucrure  is  again  quit?  efficierr.  as  the  reflected  power 
may  be  made  negligible  with  ail  the  incident  power  rrdiated  in  the  space  wre 
(refer  to  Figs.  ^.16  and  i  7).  When  £  is  greater  than  one  the  structure  upper’s 
a  surface  wave  and  again  the  power  is  transferred  from  the  space  wave  *o  the 
surface  wave  as  seen  in  Fig.  18.  Only  one  surface  wave  is  supported  witu  the 
values  of  parameters  that  permit  only  the  lowest  order  mode  to  propagate  in 
the  waveguide.  The  effect  of  varying  £,  is  shown  m  Fig.  19. 

A  check  on  the  algebra  and  computer  results  .s  possible  bv  using  the 
conservation  of  energy  principle.  The  surr.  of  the  power  radiated  in  ’be  space 
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wave,  the  power  reflected  in  the  waveguide,  and  the  power  carried  by  the  sur¬ 
face  waves  (if  any)  must  equal  the  incident  power.  This  equality  was  obtained 
to  within  0.  5  per  cent.  This  also  gave  a  check  on  the  far  field  patterns  as  the 
radiated  power  is  equal  to  a  constant  times  the  integral  of  the  square  of  the 
far  field  function.  A  verification  of  the  numerical  work  for  the  radiated  power 
is  then  an  indirect  verification  of  the  far  field  pattern. 

A  comparison  of  the  numerical  results  for  the  TEM  excitation  of  the 
surface  wave  structure,  given  in  Fig.  18,  with  those  published  by  Angulo  and 
Chang  [1959]  shows  that  they  are  not  in  agreement.  The  results  given  here 
should  compare  with  their  results  with  "h"  set  to  zero.  Their  results  show 
that  the  reflected  power  has  a  maximum  at  k^b  approximately  1.  25  with  a 
corresponding  minimum  in  the  power  radiated.  Our  results  do  not  display 
these  phenomena.  Another  paper  by  Angulo  and  Chang  [1958]  gives  the  results 
for  a  cylindric.l  geometry.  The  results  published  there  have  the  functional 
f<  -m  of  our  results  given  in  Fig.  18.  One  would  not  expect  the  change  in  the 
geometry  f’ im  cylindrical  to  rectangular  to  cause  the  change  in  response  as 
found  in  their  two  papers. 


POWER 


Fig.  17  Power  distribution  and  far  field  patters:''  i<>r  a  TKV  e \c  :  ?  ed  par- 
alle!  plate  waveguide  radiating  in  free  ••parr. 
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POWER 


16  Power  distribution  and  far  Si.'ld  patterns  for  a  TE«  excited 
surface  wave  structure. 
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Fig.  19  Power  distribution  and  far  field  patterns  for  a  TEM  excited 
structure  as  a  function  of  K. 
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5.  CONCLUSION 

The  solutions  of  three  Wiener-Hopf  type  boundary  value  problems 
have  been  given  in  this  work.  The  problem  of  a  parallel  plate  waveguide  with 
one  nlate  truncated  and  radiating  in  free  space  was  solved  in  section  2.  A  ver¬ 
ification  of  the  factorization  was  obtained  by  reference  to  the  solution  of  a  semi- 
infinite  parallel  plate  waveguide  radiating  in  free  space  given  in  Noble  [1958 L 
The  function  to  be  factored  is  the  same  for  both  problems. 

The  excitation  of  a  dielectric  slab  (surface  wave  structure)  and  the  x 
citation  of  an  incompressible,  isotropic,  plasma  slab  by  means  of  a  truncated 
parallel  plate  waveguide  were  given  in  section  3.  Both  TE  and  TEM  polariza¬ 
tions  of  the  exciting  field  were  considered.  The  results  for  the  TEM  excitation 
of  the  surface  wave  structure  were  compared  with  those  obtained  by  Angulo 
and  Chang  [1959]  and  the  differences  noted.  They  used  a  formal  factorization 
procedure  in  their  paper.  The  graphical  results  for  the  TE  excitation  of  the 
surface  wave  structure  and  the  graphical  results  for  both  TE  and  TEM  excita¬ 
tions  of  the  incompressible,  isotropic,  plasma  slab  presented  in  section  4  have 
not,  to  the  best  of  the  author's  knowledge,  been  given  elsewhere. 

The  factorization,  one  of  the  key  steps,  was  obtained  by  a  technique 
described  in  this  work.  The  factorization  was  obtained  iby'  taking  the  limit,  as 
the  transverse  dimension  approaches  infinity,  of  the  function  and  factorization 
appropriate  to  the  related  closed- region  structure.  A  closed  form  of  the  fac¬ 
torization  was  obtained  only  in  section  2.  In  the  more  difficult  problems  dis- 
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cussed  in  section  3,  the  factorization  was  in  a  form  convenient  for  numerical 
processing. 

This  technique  for  obtaining  the  factorization  is  certainly  applicable 
to  other  open- region  problems  as  discussed  in  section  2.  5.  For  example 
this  technique  should  prove  useful  for  finding  the  electromagnetic  fields  asso¬ 
ciated  with  an  incompressible,  anisotropic,  plasma  slab  when  excited  by  a 
truncated  parallel  plate  waveguide. 
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APPENDIX 


Region  of  Analyticity  of  a  Function  Given  by  ar.  Integral  Representation 
Let  G((X  )  be  defined  by  equation  (A-l) 

G(*)  =  £  J  (a.,1)  oY  (A-l) 

The  conditions  under  which  Gr(u.)  is  analytic  are  given  in  Noble  [1958],  page 
11,  and  are  presented  here  for  convenience. 

Theorem:  Let  g(  ex  ,  f  )  -  f(  ~f  )h(  a  ,  )  sat.sfv  the  conditions 

(i)  h(  a  ,  y  )  is  a  continuous  function  of  the  complex  variables  CL  and  ~J 
where  OL  lies  inside  a  region  R  and  Y  lies  on  a  contour  C. 

(li)  h(  CL  ,  Y  ^  is  a  regular  function  of  a  in  R  for  every  Y  011 

(iii)  f('J)  has  only  a  finite  number  of  finite  discontinuities  or.  C  and  a 
finite  number  of  maxima  ana  minima  on  anv  finite  part  of  C. 

(iv)  f(cf )  is  bounded  except  at  a  finite  number  of  points,  If  q  1s 

such  a  point,  so  that  ^(cL  ,  J  )  -*-  oe  as  >  ~y  then 


(  J7  -  li«  j  ’■? 

c  6 C-6  J 


exists  wheTe  the  notation  (C  -  6  )  denotes  the  concur  C  apa^t  from  a  small 
length  sum  jading  y  .  and  lim  (  6  — ►  r  i  denotes  the  limit  as  this  ex¬ 

cluded  length  tends  to  zero.  The  limit  must  be  approached  uniformly  when  CL 
lies  in  ary  closed  domain  R  within  R. 

(v)  If  C  gets  to  infinity  then  any  bounded  part  of  C  must  be  smooth 


and  renditions  (i)  and  { i : )  must  be  satisfied  for  anv  bounded  part  of  C.  The  m- 
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*e  integral  defining  G(a.)  must  be  uniformly  convergent  when  <3-  li^-s  m 
closed  domain  R‘  within  R. 

Then  G(cl)  defined  by  (A - 1 )  is  a  regular  function  of  CL  in  R. 
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13  ABSTT*CT  (continued) 

of  tne  function  and  factorization  appropriate  to  the  closed-region  structure  By 
this  means  the  factorization  and  hence  the  solut'on  to  the  open-region  boundary 
value  j.roblen  is  o^tai.  od. 

It  is  also  found  that  the  limiting  procedure  may  be  used  to  obtain  more  than 
just  the  open-region  factorization  It  is  shown  that  the  limit  of  the  complete 
closed- region  solution  becomes  the  open-region  solution.  Hence,  U.  s  yeilti:  one 
possible  method  for  the  solution  of  p  oblems  of  this  type 

The  results  of  he  numerical  imputations  are  presentee  Th  se  include  the 
average  power  reflected  in  the  waveguide,  the  average  power  radiated  in  the  space 
wave,  the  av  rage  power  transmitted  by  the  s  rface  waves,  and  tre  radiation  pa  tern 
of  the  sprje  wave 


